AD-A154  326  SOIL-AIR  INTERACTIONS  DURING  AIRBLAST-INDUCED  GROUND 
HOTIONS<U)  CALIFORNIA  RESEARCH  AND  TECHNOLOGV  INC 
CHATSWORTH  M  ROSENBLATT  ET  AL.  81  JAN  82  CRT-3470F 
UNCLASSIFIED  DNA-TR-81-103  DNA001-81-C-0247  F/G  8/7 


OTIC  FILE  COPY  AD-A154  326 


rV 


DNA-TR-81-103 


i  SOIL-AIR  INTERACTIONS  DURING 

*  AIRBLAST-INDUCED  GROUND  MOTIONS 

w1 

Martin  Rosenblatt 
Paul  J.  Hassig 
Dennis  L.  Orphal 

California  Research  &  Technology,  Inc. 

20943  Devonshire  Street 
Chatsworth,  CA  91311-2376 


1  January  1982 


Technical  Report 


CONTRACT  No.  DNA  001-81-C-0247 


APPROVED  POR  PUBLIC  RELEASE; 
DISTRIBUTION  UNLIMITED. 


THIS  WORK  WAS  SPONSORED  BY  THE  DEFENSE  NUCLEAR  AGENCY 
UNDER  RDT&E  RMSS  CODE  B34408I466  Y990AXSG00006  H2590D. 


Prepared  for 
Director 

DEFENSE  NUCLEAR  AGENCY 
Washington,  DC  20305 


DT1C 

ELECTE!^ 
MAY  3  1  1985  1 

B 


85  03 


29  009 


Destroy  this  report  when  it  is  no  longer 
needed.  Do  not  return  to  sender. 

PLEASE  NOTIFY  THE  DEFENSE  NUCLEAR  AGENCY, 
ATTN:  STTI ,  WASHINGTON,  D.C.  20305,  IF 
YOUR  ADDRESS  IS  INCORRECT,  IF  YOU  WISH  TO 
BE  DELETED  FROM  THE  DISTRIBUTION  LIST,  OR 
IF  THE  ADDRESSEE  IS  NO  LONGER  EMPLOYED  BY 
YOUR  ORGANIZATION. 


UNCLASSIFIED 


SECURITY  CLASSIFICATION  Of  THIS  PAGE 


la  REPORT  SECURITY  CLASSIFICATION 

UNCLASSIFIED 


2 a  SECURITY  CLASSIFICATION  AUTHORITY 

N/A 


2b  DECLASSIFICATION/  DOWNGRADING  SCHEDULE 

N/A  since  UNCLASSIFIED 


a  PERFORMING  organization  REPORT  NUMBER(S) 

CRT3470F 


6a  NAME  OF  PERFORMING  ORGANIZATION  6b.  OFFICE  SYMBOL 

California  Research  &  of  applicable) 

Technology,  Inc. 


6c.  ADORESS  (City,  Stare,  and  ZIP  Code) 

20943  Devonshire  Street 
Chatsworth,  California  91311-2376 


REPORT  DOCUMENTATION  PAGE 


lb.  RESTRICTIVE  MARKINGS 


3  DISTRIBUTION /AVAILABILITY  OF  REPORT 

Approved  for  public  release;  distribution 
is  unlimited. 


S.  MONITORING  ORGANIZATION  REPORT  NUMBER(S) 

DNA-TR-81-103 


7a.  NAME  OF  MONITORING  ORGANIZATION 

Director,  Defense  Nuclear  Agency 


7b.  ADDRESS  (City,  State,  and  ZIP  Code) 

Washington,  DC  20305 


8a.  NAME  OF  FUNDING/SPONSORING 
ORGANIZATION 


8b  OFFICE  SYMBOL 
(If  applicable) 


9  PROCUREMENT  INSTRUMENT  IDENTIFICATION  NUMBER 

DNA  001-81-C-0247 


c.  ADDRESS  (G'fy,  Sfafe,  and  ZIP  Code) 

10.  SOURCE  OF  FUNDING  NUMBERS 

PROGRAM 
ELEMENT  NO 

62715H 

PROJECT 

NO 

Y990AXS 

TASK 

NO 

G 

11  title  (Include  Security  Classification ) 

SOIL-AIR  INTERACTIONS  DURING  AIRBLAST-INDUCED  GROUND  MOTIONS 


12  personal  author(S)  Martin  Rosenblatt 
Paul  Joseph  Hassig 


13a  TYPE  OF  REPORT 

Technical  Report 


nblatt  Dennis  Lee  Orphal 

Hassig 


13b  TIME  COVERED  14  OATS  Of  REPORT  (Year,  Month,  Day)  Il5  PAGE  COUNT 

FROM  81  Jul  TO  82  Jan  |  1982,  January  ]  160 


'6  SUPPLEMENTARY  NOTATION 

This  work  sponsored  by  the  Defense  Nuclear  Agency  uncier  RDT&E  RMSS 
Code  B344081466  Y99QAXSG00006  H2590D. 


17  COSATi  CODES _ I  18  SUBJECT  TERMS  ( Continue  on  reverse  if  necessary  and  identify  by  block  number) 

su8-group  |  Ground  Motion  Multiphase  Flow  MISERS  BLUFF 

Airblast  Two-Phase  Flow 

_  Pore-Air Fluidization 


19  abstract  (Cont»nue  on  reverse  if  necessary  and  identify  by  block  number) 

Some  ground  motion  records  in  HE  field  tests,  e.g.,  MISERS  BLUFF  1 1 -2  and  PRE-MINE 
THROW  IV-6,  exhibit  an  enhanced  upward  velocity  immediately  following  the  initial 
airblast-induced  downward  motion.  Also,  laboratory  experiments  show  the  development  of 
upward  soil  motions  during  positive  air  overpressure  loadings;  the  motion  in  the 
laboratory  tests  was  shown  to  be  related  to  the  flow  of  the  air  through  the  soil  pores. 

In  this  report,  the  effects  of  pore-air*  on  ground  motions  in  dry  porous  soils  are 
analyzed  using  five  numerical  simulations  in  one  space  dimension  and  time.  Also, 

*Dry  porous  soils  can  have  30-40%  of  their  volume  filled  with  air;  this  in-situ  air  is 
called  pore-air  and  is  one  important  aspect  of  soil/airblast  interactions. 


20  DISTRIBUTION /AVAILABILITY  OF  ABSTRACT 

□  UNCLASSIFIED/UNLIMITED  C3  SAME  AS  RPT.  □  OTIC  USERS 


22a  NAME  OF  RESPONSIBLE  INOIVIOUAL 

Betty  L.  Fox 


OO  FORM  1473,  84  MAR  83  APR  edition  may  be  used  until  exhausted 

All  other  editions  are  obsolete 


21  ABSTRACT  SECURITY  CLASSIFICATION 

Unclassified 


22b.  TELEPHONE  (Include  Area  Code)  22c  OFFICE  SYMBOL 

202-325-7042  DNA/STTI 


SECURITY  CLASSIFICATION  OF  THIS  PAGE 


UNCLASSIFIED 


ir  ^  .  J" 


ST.  VJ  f 


UNCLASSIFIED 

SECURITY  CLASSIFICATION  OF  THIS  PAGE 


19.  ABSTRACT  (Continued) 

experimentally  observed  ground  motions  involving  enhanced  vertical  velocities  are 
explained  using  two-phase  (soil  and  air)  physics,  including  (a)  the  flow  of  air 
through  the  soil  pores,  and  (b)  soil  fluidization.  Fluidization  of  the  soil  occurs 
when  the  soil  particles  lose  contact  with  each  other  and  the  total  pressure  is 
entirely  supported  by  the  pore-air. 

The  primary  conclusions  are: 

1.  Large  sustained  upward  velocities  can  develop  in  near-surface  soils 
subjected  to  airblast  loading  due  to  soil  fluidization  caused  by  pore-air 
phenomena. 

2.  In  the  laboratory  experiments  involving  relatively  high  soil  permeabilities, 
the  100  psi  peak  overpressure  airblast  caused  air  to  flow  through  the  soil 
pores,  thereby  increasing  pore-air  pressures  which  resulted  in  soil 
fluidization  during  the  positive  overpressure  phase.  Peak  pore-air  flow 
speeds  of  -15  m/s  are  predicted,  and  the  soil  "fluidization  front" 
propagates  at  ~4  m/s  (Case  1). 
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in  the  soil  expanded  to  reach  pressure  equilibrium  with  the  low  airblast 
pressure. 

4.  The  DICE  code  numerical  simulations  (Cases  1  to  5)  reproduced  the 
phenomenology  observed  in  the  laboratory  and  field  experiments.  The 
relatively  simple  models  used  for  soil  permeability  and  internal  interaction 
forces  (drag)  were  adequate  for  quantitatively  reproducing  the  experimental 
data;  however,  the  iri-situ  soil  permeability  is  a  sensitive  physical 
parameter  and  careful  measurements  and  model  calibrations  will  be  necessary 
before  quantitative  predictions  and  associated  uncertainties  of  ground 
motions  can  be  made  concerning  strategic/tactical  sites  of  interest. 
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SECTION  1 


INTRODUCTION  AND  SUMMARY 


1.1  BACKGROUND 

The  airblast  from  near-surface  explosions  can  induce  motions 
in  the  ground  which  will  damage  structures  near  the  surface. 
Predictive  calculations  of  airblast-induced  ground  motion  have 
generally  assumed  that:  (1)  the  air/ground  interface  acts  as  an 
impermeable  membrane  and  (2)  any  effects  of  air  contained  in  the 
ground  are  incorporated  in  the  pore  crush-up  portions  of  the 
constitutive  equations  used  to  model  the  stress-strain  response 
of  the  ground  material. 

Observations  and  measurements  from  several  field  and 
laboratory  experiments  suggest  that  a  more  complex  multiphase 
interaction  between  the  air  and  soil  occurs  which  may  substan¬ 
tially  influence  near-surface  ground  motion  wave  forms.  In 
particular,  some  vertical  velocity  records  from  several  HE  field 
tests,  e.g.,  MISERS  BLUFF  II-21  and  PRE-MINE  THROW  IV-62,  exhibit 
an  enhanced  upward  velocity  immediately  following  the  initial 
airblast-induced  downward  motion.  Also,  laboratory  experiments 
show  the  development  of  upward  soil  motions  during  positive  air 

3 

overpressure  loading  ;  this  motion  was  shown  to  be  related  to  the 
flow  of  the  air  through  the  soil  pores. 

Porosities  of  soils  may  be  30-40%,  which  means  that  the 
solid  soil  particles  form  a  "lattice"  or  connected  network  which 
contains  30-40%  pore  volume  filled  with  air  and/or  water.  For 
dry  soils,  pore-air  effects  refer  to  phenomena  related  to  the 
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aoil  pores  being  filled  with  air  and  interacting  with  the  solid 
soil  particles. 

In  this  study,  the  effects  of  pore-air  on  ground  motions  in 
dry  porous  soils  are  numerically  simulated  in  one  space  dimension 
and  time.  Also,  experimentally  observed  ground  motions  involving 
enhanced  vertical  velocities  are  explained  using  two-phase  (soil 
and  air)  physics  including  the  flow  of  air  through  the  soil 
pores,  and  fluidization  of  the  soil. 

1.2  OBJECTIVES  AND  APPROACH 

4 

The  numerical  simulations  use  the  multiphase  DICE  code 
which  has  previously  been  applied  to  nuclear  burst  problems 
involving  soil-air  interactions  above  the  ground  surface.*  The 
mixing  of  solid/liquid  soil  particles  with  the  atmosphere  and 
fireball  generally  involved  dust/pebble  concentrations  which  were 
comparable  to  the  air  density  (p~10-3  gm/cm3  )  or  lower.  Thus, 
most  of  the  soil-air  flow/ interact ions  were  in  the  aerodynamic 
drag  regime.  However,  the  current  investigation  involves  the 
near-surface  soil  material  which  has  an  initial  soil  density  of 
~2  gm/cm3. 

The  research  objectives  were  to: 

1.  Develop  and  implement  a  consistent  flow/ inter act ion 
(i.e.,  generalized  drag)  model  involving  soil  and 
pore-air  constituents  for  all  soil  to  air  density 
rat ios . 


*The  so il  material  was  injected  into  the  atmosphere  as  dust/pebble 
particles  using  analytical  models  for  crater  ejecta  (in  surface 
bursts)  or  surface  sweep-up  of  soil  particles  (in  air  bursts). 
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2.  Numerically  simulate  and  analyze  the  pore-air 
phenomenology  and  induced  ground  motions  in  one  space 
dimension  at  the  100  psi  peak  overpressure  range  after  a 
1  Mt  burst  on  porous  soils  with  various  permeabilities. 

3.  Perform  and  analyze  cases  which  could  be  compared  with 
laboratory  and  field  test  experimental  data. 

A  model  describing  the  soil-air  interaction  and  flow  has 
been  developed  and  implemented  in  the  multiphase  DICE  code.  The 
flow/ inter act ion  model  specifies  the  internal  mutual  interaction 
forces  between  the  soil  "lattice"  and  the  pore-air  as  a  function 
of  relative  velocity,  soil  porosity,  and  a  characteristic  soil 
particle  size.  The  model  is  described  in  Section  2.3. 

Table  1.1  summarizes  the  five  cases  performed  and  analyzed. 
In  Cases  1  through  4,  vertical  ground  motions  are  numerically 
simulated  for  homogeneous  porous  soils  which  vary  in  initial 
permeability*  by  a  factor  of  8.  The  airblast  history  (100  psi 
peak  overpressure)  applied  to  the  soil  surface  is  similar  to  a 

5 

1  Mt  ideal  airblast  history  (see  the  bottom  of  Table  1.1).  Case 
5  is  a  calculation  of  the  vertical  ground  motions  corresponding 
to  one  location  in  the  HE  multiburst  experiment  MISERS  BLUFF  1 1-2 
(see  Section  4).  Three  ground  material  layers  were  modeled  based 
on  measured  material  properties  including  hysteretic  load-unload 
paths.  The  measured  airblast  loading  used  in  the  calculation  is 
also  shown  on  Table  1.1. 

Cases  1  to  4  are  described  and  compared  to  laboratory 

3 

vertical  shock  tube  tests  in  Section  3,  and  Case  5  is  described 

*Note  that  the  soil  permeability  (k)  depends  on  both  porosity  and 
the  characteristic  particle  size  (Dp)  as  described  in  Section 
2.3.  Therefore,  Dp  will  be  used  to  show  differences  in 
permeability  for  tne  various  cases,  where  k  a  D2  . 
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and  compared  to  the  MISERS  BLUFF  measured  ground  motions  in 
Section  4.  Multiphase  flow  theory  and  DICE  code  modeling 
techniques  are  described  in  Section  2. 


1.3  SUMMARY  OF  RESULTS 

3 

Pore-air  phenomena  were  studied  by  Zernow  et  al.  in  a 
series  of  laboratory  experiments.  The  experiments  involved  the 
use  of  a  shock  tube  to  load  the  surface  of  a  column  of  sand  with 
overpressure  histories  similar  to  the  positive  phase  of  the  ideal 
airblast  at  the  100  psi  range  from  a  1  Mt  nuclear  burst.  The 
motion  of  the  sand  column  was  observed  photographically.  The 
DICE  numerical  simulation  Cases  1  to  4  had  a  similar  overpressure 
loading  function  applied  at  the  soil  surface. 


Figure  1.1  shows  the  surface  displacements  observed  in  the 
laboratory  experiments  compared  to  those  calculated  in  DICE  Cases 
1  and  2.  The  1.78"  diameter  orifice  experiments  (see  Section  3) 
had  a  faster  decay  rate  in  the  airblast  loading  compared  to  the 
1.36"  diameter  orifice  experiments,  as  shown  below: 
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Figure  1.1.  Main  Surface  Displacement  Versus  Time  Observed  in 
Experiments  (Zernow  ef  at .,  1973)  and  Calculated 
by  DICE  (Cases  1  and  2). 


The  result  is  an  earlier  upward  motion  of  the  main  sand  column 
(-.200  s  compared  to  -.400  s)  and  higher  velocity  (-75  cm/s 
compared  to  -45  cm/s)  for  the  1.78"  diameter  orifice  experiments. 
Maximum  heights  of  nearly  20  cm  are  attained  by  the  lofted  sand. 
In  the  absence  of  a  sustained  upward  force  acting  on  the  sand, 
the  initial  lofting  velocity  needed  to  reach  20  cm  height  is 
V  =  /2gH  -200  cm/s.  Thus,  the  experimental  and  numerical  results 
provide  convincing  evidence  of  the  sustained  application  of  an 
upward  force  associated  with  pore-air  phenomena. 

The  surface  displacement  in  Case  1  is  similar  to  the  1.78" 
diameter  orifice  experiment.  Both  surfaces  start  to  move  upward 
at  t  -.200  s  and  maintain  an  approximately  constant  upward 
velocity  until  t  -.400  s.  The  continuing  rise  of  the  soil 
surface  in  Case  1  above  the  experimental  curve  is  due  to  the 
longer  duration  of  the  airblast  loading  function  used  in  the 
numerical  simulation. 

The  upward  surface  displacement  in  Case  2  also  begins  at 
t  -.200  s,  but  a  smaller  constant  upward  velocity  develops.  The 
smaller  displacements  and  velocities  compared  to  Case  1  are 
indicative  of  a  smaller  sustained  upward  force  resulting  from  the 
permeability  dependence  of  the  internal  interaction  (drag) 
forces.  The  factor  of  2  smaller  characteristic  particle  size 
used  in  Case  2  yields  a  factor  of  4  decrease  in  modeled 
permeability  compared  to  Case  1. 

The  fact  that  the  soil  surface  displaces  upward  during  the 
positive  air  overpressure  loading  on  the  surface  is  a  pore-air 
flow  phenomenon.  High  pressure  air  flows  into  the  porous  soil 
from  the  airblast  loading  on  the  surface.  As  air  flows  or 
permeates  downward  through  the  soil  pores,  the  pore-air  pressure 
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increases  and  the  total  applied  pressure  load  supported  by  the 

solid  soil  particles  decreases.  Thus,  P.  .  ,  «  P  P.. 

*  total  soil  air 

Fluidization  of  the  soil  occurs  when  the  soil  particles  lose 
contact  with  each  other  and  the  total  pressure  is  entirely 
supported  by  the  pore-air .  Once  fluidization  occurs,  the 
fluidized  soil  displaces  upward  (with  respect  to  the  unfluidized 
soil).  These  phenomena  are  described  in  more  detail  in  Section 
3.2  for  Cases  1  to  4  and  illustrated  on  Figures  1.2  and  1.3. 

Figure  1.2  shows  the  pressure  and  velocity  versus  depth 
profiles  at  t  -  .200  s  for  Cases  1  and  2.  The  pressure  profiles 
(at  top  of  figure)  show  the  extent  of  fluidization  to  be  ~40  cm 
in  Case  1  and  ~25  cm  in  Case  2.  The  velocity  profiles  (at  bottom 
of  Figure  1.2)  indicate  that  there  is  an  upward  acceleration  of 
the  soil  material  associated  with  fluidization.  The  soil  is 
moving  downward  at  a  constant  velocity  of  -  23  cm/s  below  the 
fluidized  region.  The  air  is  permeating  down  through  the  soil 
below  the  fluidized  region,  and  is  flowing  up  and  out  of  the  soil 
in  the  fluidized  region.  For  the  more  permeable  soil  modeled  in 
Case  1,  the  air  has  larger  downward  velocities  and  flows  through 
the  soil  pores  to  greater  depths.  A  greater  relative  velocity 
between  the  air  and  soil  is  also  present  in  the  fluidized  region 
for  Case  1  as  compared  to  the  less  permeable  soil  in  Case  2. 

Figure  1.3  shows  the  pore-air  overpressure  and  soil  particle 
velocity  versus  time  at  various  depths  in  the  soil  for  Case  1. 

At  a  given  depth,  the  air  overpressure  gradually  rises  to  a  peak 
and  then  decays  to  a  level  approximately  equal  to  the  applied 
surface  overpressure  at  which  time  the  soil  becomes  fluidized. 

The  soil  velocity  histories  (Figure  1.3b)  show  that  up  to 
t  -.100  s  the  entire  soil  column  essentially  moves  downward  as  a 
rigid  body.  Fluidization  at  each  station  depth  becomes  apparent 
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Figure  1.2.  Pressure  and  Velocity  Profiles  Versus  Depth  at  t  =  .200  Seconds  for 
Cases  1  and  2. 


after  this  time  when  the  soil  flowing  at  that  depth  suddenly 
accelerates  its  velocity  above  the  velocity  of  the  soil  column 
below  that  depth. 

The  progression  of  fluidization  down  through  the  soil  column 
is  shown  on  Figure  1.4a  for  Cases  1  to  4.  For  Case  1  the 
"fluidization  front"  moves  down  the  soil  at  ~4  m/s.  The 
fluidization  front  moves  more  slowly  at  ~2  m/s  for  the  less 
permeable  soil  modeled  in  Case  2.  In  Case  3,  the  soil  at  a  given 
depth  is  fluidized  earlier  than  in  Case  2  because  there  is  no 
lithostatic  pressure  distribution  (g  =  0)  in  the  soil  to  be 
overcome  by  the  pore-air  before  fluidization  occurs.  A  fluidiza¬ 
tion  front  velocity  of  ~7  m/s  is  predicted  for  the  highly 
permeable  soil  in  Case  4. 

Figure  1.4b  shows  the  calculated  peak  air  overpressure 
arrival  times  versus  depth  for  Cases  1  to  4;  the  experimental 
data  are  also  shown  for  comparison.  The  more  permeable  soils 
allow  air  to  travel  faster  down  the  soil  pores,  causing  earlier 
arrival  times  for  the  peak  ovepressure.  The  slope  of  the  peak 
arrival  time  versus  depth  curves  for  Cases  1  to  4  are  quite 
similar  to  the  soil  fluidization  time  versus  depth  curves  in 
Figure  1.4a.  For  each  case  there  is  a  delay  of  ~.080  to  .120  s 
after  the  arrival  of  the  peak  overpressure  at  a  given  depth  until 
soil  fluidization  begins  at  that  depth. 

The  MISERS  BLUFF  1 1-2  experiment  involved  the  simultaneous 
detonation  of  six  high-explosive  charges  forming  a  hexagon  100  m 
on  a  side.  An  approximation  to  the  airblast  measured  at  25  m 
radius  from  the  center  of  the  hexagon  was  used  in  Case  5  as  the 
overpressure  loading  function.  In  this  case,  the  airblast 
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included  a  negative  overpressure  phase  beginning  at  t  -.150  s  and 
returning  to  zero  overpressure  at  t  -.400  s. 

The  material  properties  for  the  MISERS  BLUFF  1 1-2  calcula¬ 
tion  (Case  5)  involved  three  soil  layers  including  hysteresis  and 
the  associated  variation  in  load-unload  constrained  moduli.  The 
initial  soil  permeabilities  in  Case  5  were  much  smaller  than  in 
Cases  1  to  4  where  the  initial  permeabilities  were  in  the  10  2  - 
103  Darcy  range  (1  Darcy  -  9.87  x  10~9  cm2).  In  Case  5,  the 
initial  permeabi it ies  were  different  in  each  soil  layer  -  k  =  .04 
Darcy  for  the  top  layer,  k  *  69  Darcy  for  Layer  2  and  k  -  .003 
Darcy  for  Layer  3.  The  low  permeability  in  the  top  layer  limited 
the  amount  of  high  pressure  air  which  flowed  into  the  soil. 

Also,  the  internal  interaction  forces  between  the  air  and  soil 
particles  (i.e.,  drag)  increased  with  the  decrease  in 
permeability,  so  that  the  relative  air-soil  velocities  were 
nearly  zero. 

Sustained  upward  velocities  and  fluidization  of  the  MISERS 
BLUFF  1 1 -2  soils  occurred  primarily  because  of  the  negative 
overpressures  in  the  air blast  and  the  associated  expans  ion  of  the 
pore-air  in  the  soil;  the  pore-air  expands  so  that  pressure 
equilibrium  can  be  attained  with  the  applied  airblast  pressure. 

If  the  soils  were  not  hysteretic,  this  fluidization  would  begin 
when  the  negative  overpressures  were  first  applied  to  the 
incohesive  soils.  The  hysteresis  in  the  soils,  however,  permit 
fluidization  to  occur  at  small  positive  overpressures  since  the 
soil  irreversibly  compacts  and  cannot  support  stresses,  upon 
unloading,  at  soil  densities  slightly  above  the  initial  density. 

Various  ground  motion  gauges  originally  at  0.5  m  depth  and 
2b  m  radius  in  MISERS  BLUFF  1 1-2  measured  the  vertical  velocity 
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(and  thus  displacement)  time  histories.  Figure  1.5  compares  the 
experimental  vertical  velocity  and  displacement  data  with  the 
Case  5  predictions  of  a  tracer  particle  originally  at  0.5  m 
depth.  The  sustained  upward  velocity  from  t  -.150  s  to  beyond 
t  -.400  s  demonstrates  the  important  role  of  pore-air  phenomena 
during  the  negative  overpressure  phase  in  low  permeability  soils. 

1.4  CONCLUSIONS  AND  RECOMMENDATIONS 

1.  Large  sustained  upward  velocities  can  develop  in 
near-surface  soils  subjected  to  airblast  loading  due  to  soil 
fluidization  caused  by  pore-air  phenomena. 

3 

2.  In  the  laboratory  experiments  involving  relatively  high 
soil  permeabilities,  the  100  psi  peak  overpressure  airblast 
caused  air  to  flow  through  the  soil  pores,  thereby  increasing 
pore-air  pressures  which  resulted  in  soil  fluidization  during  the 
positive  overpressure  phase.  Peak  pore-air  flow  speeds  of 

-15  m/s  are  predicted,  and  the  soil  fluidization  front  propagates 
at  -4  m/s  (Case  1) . 

3.  In  the  MISERS  BLUFF  1 1 -2  field  experiment1  involving  low 
soil  permeabilities  and  almost  no  pore-air  flow,  soil  fluidiza¬ 
tion  and  sustained  upward  velocities  occurred  during  the  negative 
overpressure  phase  as  the  pore-air  in  the  soil  expanded  to  reach 
pressure  equilibrium  with  the  low  airblast  pressure. 

4.  The  DICE  code  numerical  simulations  (Cases  1  to  5) 
reproduced  the  phenomenology  observed  in  the  laboratory  and  field 
experiments.  The  relatively  simple  models  used  for  soil 
permeability  and  internal  interaction  forces  (drag)  were  adequate 
for  quantitatively  reproducing  the  experimental  data;  however. 
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Figure  1.5.  (a)  Vertical  Velocity  and  (b)  Vertical  Displacement 

Time  History  Comparisons,  Case  5  Versus  MISERS  BLUFF  I 1-2 
Experimental  Data  at  25  m  Radius  and  0.5  m  Initial  Depth. 


the  in-situ  soil  permeability  is  a  sensitive  physical  parameter 
and  careful  measurements  and  model  calibrations  will  be  necessary 
before  quantitative  predictions  and  associated  uncertainties  of 
ground  motions  can  be  made  concerning  strategic/tactical  sites  of 
interest. 


5.  The  following  recommendations  are  suggested: 


•  Evaluate  the  pore-air  and  associated  soil  fluidization 
effects  on  ground  motions  for  peak  overpressures  up  to 
1  kbar  ~ 15, 000  psi. 

•  Obtain  in-situ  soil  permeability  measurements  for 
strategic  and  tactical  sites  of  interest. 

•  Extend  the  pore-air  analysis  capability  to  pore-air -water 
environments  since  most  soil  pores  contain  some  water; 
and  in  fact,  many  cases  of  interest  involve  saturated  or 
nearly  saturated  soils  below  some  depth. 

•  Analyze  the  importance  of  pore-air  and  soil  fluidization 
phenomena  in  two  space  (2-D)  dimensions  for  propagating 
airblast  waves  which  interact  with  surface  dust/pebble 
material.  (See  following  sketch.)  Both  experimental  and 
theoretical  analyses  should  be  conducted. 


Schematic  of  Two-Dimensional  Soil-Air 
Interactions  Behind  Air  Shock 
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SECTION  2 


MULTIPHASE  THEORY 


2.1  GOVERNING  PHYSICAL  EQUATIONS  -  GENERAL  FORMULATION 


Consider  a  bounded  volume  fl  with  surface  area  S  which 
contains  physically  distinct  material  phases.  Distinct  phases 
can  include  solid,  liquid  and/or  vapor  phases  as  well  as 
individual  phases  distinguished  by  particle  type/size.  If  the 
volume  of  each  phase  is  n^,  where  k  =  1  to  K,  then 


n  =  E  n 

k=l,K  K 


or 


1  -  E  n.  /n  =  E  » 

k=l,K  k=l,K 


where 


v  «  n,  /n  is  the  Volume  Fraction. 

J*k. 


(2.1) 


The  continuum  mechanics  concept  of  a  mass  dens itv  for  each  phase 
is  assumed.  Thus,  as  a  test  volume  n  and  associated  phase 
volume  become  small,  the  mass  Mk  and  the  average  and  local 
densities  are  related  by 

Pk  -  lim  (MZ/fO  ;  pR  is  average  density 

n  — o 

Pk  -  lim  (M^/nk)  -  Pk/vk  ?  Pk  is  local  density 


The  total  mass  density  is  then 


total 


k*l,K 


E  i 

k=l,K 


pkvk 


(2.2) 


The  dynamic  physical  equations  for  each  phase  (and  thus  for 
the  entire  multiphase  system)  involve  identical  considerations 
concerning  conservation  of  mass,  momenta  and  energy.  For 
simplicity,  and  without  any  loss  in  generality,  the  subscript  k 
for  each  phase  will  be  dropped  in  the  following  derivation. 

Thus,  the  mass,  momenta  and  energy  of  any  one  material  phase  in  a 
volume  n  are  given  by 


M(t)  =  /pdV 
n 


MOM  (t) 


fp u  dV  *  /pudV 

n  n 


(2.3) 


EN(t )  -  /p(e  +  ^u^u  )dV 

n  ] 


/pEdV 

a 


where 


p  is  the  average  dens ity  of  one  phase  (p  =  pv), 
u1  is  the  local  velocity  vector  of.  one  phase . 
e  is  local  specif ic  internal  energy  of  one  phase,  and 
E  is  specif  ic  total  energy  ■  e  +  1su1u^  of  one  phase. 


The  summation  convention  on  repeated  indices  is  assumed,  and  a 
comma  preceding  an  index  will  indicate  covariant  differentiation. 
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The  dynamic  equations  of  motion  will  be  derived  for  arbi¬ 
trary  changes  in  the  volume  Q.  -  Q(t)  since  the  DICE  code  uses  an 
expanding  computational  grid  which  is  not  coupled  to  the  material 
particle  velocity-  The  horizontal  and  vertical  grid  lines  are 
permitted  to  move  with  arbitrary  velocities,  Ug .  Thus,  in 
general,  the  volume  n  and  surface  area  S  of  the  computational 
cells  will  vary  with  time  in  a  prescribed  fashion.  The  following 
fundamental  equation  relates  the  time  derivative  of  a  volume 
integral  to  the  arbitrary  surface  velocity  Ug  (Ref.  6,  Eqn.  81.4) 

d/dt  St dV  =  fdf/Zt  dV  +  jyu^dA.  (2.4) 

n  n  s  b  3 

where  f  =  iTCx^t)  is  an  arbitrary  tensor  and  dA_.  is  the  outward 
surface  vector. 

Equations  (2.3)  and  (2.4)  can  be  used  to  obtain  the  dynamic 
mass,  momentum  and  energy  equations  in  terms  of  volume  and 
surface  integrals.  The  approach  for  the  mass  equation  is 
indicated  below. 

Letting  f  *  p  in  Equation  (2.4)  gives 

dM(t)/dt  =  fdp/dt  dV  +  /pu3dA.  (2.5) 

n  s  b  3 

The  usual  global  conservation  of  mass  equation  in  integral 
form  is  obtained  by  letting  uA  -=  uA,  i.e.,  the  surface  S  moves 
with  the  particle  velocity  u  and  therefore  the  mass  contained  in 
n(t)  is  constant  if  there  are  no  mass  sources.  Thus, 

0  =  fdp/dt  dV  +  /pu3dA . 

n  s  3 


(no  mass  sources  in 
Lagrangian  frame) 


If  there  are  mass  sources  for  this  phase,  for  example  from  the 
vaporization  of  a  solid  phase  to  <.  apor  phase,  then  the  mass 
source  term  can  be  written  as 


dM(t)/dt  =  /p+dV  (mass  source  term  in  (2.6a) 

fl  Lagrangian  frame) 

where  p+  is  the  mass  source  per  unit  volume.  Thus,  one  form*  of 
the  dynamic  mass  equation  is 

/p  dV  =  /3p/3t  dV  +  /pu^dA.  (Lagrangian  frame)  (2.6b) 

n  n  s  J 


Subtracting  Equation  (2.6b)  from  Equation  (2.5)  and  rearranging 
terms,  we  obtain 

dM(t)/dt  5  d/dt  /pdV  -  /p(uj  -  uj)dA.  +  / p  dV  (2.7) 

n(t)  S(t)a  J  n(t) 

Equation  (2.7)  is  the  integral  form  of  the  mass  equation  for 
an  arbitrary  moving  volume  fi(t). 

In  a  similar  fashion,  the  momenta  and  energy  equations  are 
obtained.  However,  surface  and  body  forces  acting  on  S(t)  and  in 
fl(t)  must  now  be  considered  since  these  forces  generate  momenta 
and  energy.  The  surface  forces  per  unit  area  are  obtained  from  a 
material  stress  tensor  o^ ,  which  for  hydrodynamic  problems  is 


*Using  the  Gauss  divergence  theorem  to  change  the  surface 
integral  to  a  volume  integral,  we  find  the  usual  differential 
equation  forms 

p+  -  3p/3t  +  (pu^),^  -  p  +  pu^ 
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H  M  1 

a  J  -*  -PS  J ;  thus,  the  surface  force  Fg  is 


=  /cx^dA^ 


Surface  Force 


(2.8) 


Note  that  only  the  stresses  for  the  one  phase  being  considered 
are  integrated  to  get  F^  for  that  phase.  See  Section  2.2  for 
further  discussions  of  the  partitioning  of  stresses. 


Body  forces  are  of  two  types: 


•  External  Gravitational  Body  Force 
/pgAdV 

n 

•  Internal  Body  Force  from  the  relative  motions  and 
interactions  with  other  phases  (-drag  force) 

n  int 


Thus,  the  body  force  F^  is 

13 


rB  -  ^  *  £int>dV 


Body  Force 


(2.9) 


The  surface  and  body  forces  in  Equations  (2.8)  and  (2.9) 
represent  momenta  source  terms;  in  addition,  there  may  be 
additional  momenta  sources  due  to  phase  changes  [analogous  to  the 
mass  sources  in  Equation  (2.6a)].  Thus,  in  the  Lagrangian  or 
material  frame  of  reference,  i.e.,  u£  =  u1,  we  have 

dMOM1/dt  -  f£  +  F^  +  Jp  uidV 
S  8  n  +  + 


J 


■1  »  *  * 
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where  p+ u^  are  the  momenta/volume  source  terms  due  to  phase  changes 
And  finally,  the  integral  form  of  the  momenta  equation  for  an 
arbitrarily  moving  volume  fl(t)  is 


dMOMX/dt  =  d/dt  /pu1dV  =  /{puA(u^  - 


n(t) 


S(t) 


u^  ) 


+  o 


i} 


)dA. 


(2.10) 


fipg1  +  f*  +  p  uJdV 

n(t)  c 


Similarly,  the  energy  equation  is 


dEN/dt  =  d/dt  /pEdV  =  /{pE(u|  -  u3)  +  o^u  .  }dA. 

n(t)  S(t)  1  ] 


+  /  (pg1u .  +  f^  fcu.  +  Q+)dV 

n(L)  1  inc  1 


(2.11) 


where  Q+  represents  the  enexgy/voluroe  source  terms  due,  in 
general,  to  any  combination  of  phase  changes,  heat  conduction  and 
thermal  radiation  transport. 


Equations  (2.7),  (2.10)  and  (2.11)  represent  the  general 
dynamic  equations  governing  each  phase  in  a  mixed  phase  flow 
configuration.  Before  a  solution  to  a  specific  problem  can  be 
obtained,  the  initial  and  boundary  conditions  must  be  specified 
as  well  as  a  complete  set  of  constitutive  equations  including: 


•  Material  Equations  of  State  (elastic-plastic-hydro 
description  of  o1-*) 


Internal  Interaction  Forces  (f1  .)  as  a  function  of 
relative  velocity,  pressure,  pirosity,  etc. 
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The  generalized  approach  for  material  equations  of  state  in 

two-phase  media  is  described  in  Section  2.2,  and  the  specific 

material  parameters  are  described  in  Sections  3  and  4  with  th 

case  descriptions  and  calculat ional  results.  The  formulation 

used  for  f f  .  for  all  the  cases  is  described  in  Section  2.3. 
int 

And,  Sections  2.4  and  2.5  contain  DICE  test  cases. 

2.2  STRESS  PARTITIONING  IN  A  SOIL-AIR  MEDIA 


A  small  portion  of  a  two-phase  material  consisting  of  a 
solid  soil  phase  and  a  pore-air  phase  is  schematically  shown 
below: 


If  a  loading  force,  F^oa(j,  i®  applied  to  a  surface,  S,  then 


load 


load  1 


(2.12 


The  partitioning  of  the  load  stress,  °ioad'  soil  and 

pore-air  parts  is  accomplished  by  introducing  the  following 
equation  and  associated  variables: 


ii  -  cii 

O  J  =  0  Jy  —  p  Q  J 1/ 

load  S3  a  a 


(2.13) 


where 


v  and  v  are  the  soil  and  air  area  fractions  (V  +  v  =  1) ; 
s  a  - -  s  a 


these  are  often  assumed  to  be  equal  to  the  volume  fractions,  but 
anisotropic  pore/channeling  configurations  are  possible  which 
would  invalidate  this  assumption. 


and  61-*  are  the  local  stresses*  such  that 
s  a  - - 


/31]v  dA,  is  the  force  contribution  from  the  soil,  and 


s  s  j 


A  M 

i-p  5  J v  dA .  is  the  force  contribution  from  the  air. 
a  a  j 


S 


The  mater ial  equations  of  state  (EOS)  for  the  air  and  the 
soil  complete  the  specification  of  Equation  (2.13).  The  air  EOS 
is  relatively  straightforward: 


Pa  -  W  ea> 


General  EOS 


(2.14) 


*These  local  stresses  6^  (force/area  for  one  phase)  are  analogous 
to  the  local  density  (mass/volume  for  one  phase).  Note  that 
the  particle  velocity  and  specific  internal  energy  are  always 
local  variables  in  the  current  formulation  and  a  u  or  e  is 
suppressed . 


For  example, 


Fa  '  (r  '  1)paea 


Ideal  Gas  EOS 


The  solid  soil  stress,  Og-1,  is  generally  not  measured;  however, 
the  measured  soil  stress  can  be  identified  by  noting  that  for 
relatively  incompressible  soil  particles,  the  pore-air  pressure 
(P  )  and  the  stresses  supported  by  the  touching  so_il  particles 

gL  . 

(o^3)  both  contribute  to  the  local  soil  stress, 


a 


13  _ 


P  6 
a 


ij 


(2.15) 


and  plugging  this  equation  into  Equation  (2.13),  we  find 


;  ^  ^  V  —  £  0^"^ 

load  s  s  a° 


(2.16) 
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Thus  experiments  measuring  °^oa(j  performed  on  soils  with  pore 

ide  the  soil  EOS  (oi3),  where 


pressure  P&  will  provide 


°i]  =  and 


ot3  „  .  a13  -  P  O13 
load  s  a 


(2.17) 


The  measured  soi_l  EOS  (oi3)  can  be  a  complicated  function  of  many 
variables,  including  dem 
strain,  and  strain  rate. 


variables,  including  density  (p  ),  specific  internal  energy  (e  ) , 

S  8 


2.3  INTERNAL  INTERACTION  (OR  DRAG)  FORCE 

The  internal  body  force  per  unit  volume,  f*  t,  for  the 
simple  two-phase  system  consisting  of  air  and  soil  particles  of 
radius  R  is  modeled  in  this  subsection.  Note  that  in  a  two-phase 


\.  --V 


-1 


■;-:v 


yy  yl 

J#_. 
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system,  the  internal  force  acting  on  the  soil  by  the  air  is 
exactly  opposite  to  the  internal  force  acting  on  the  air  by  the 
soil . *  Thus , 


(2.18) 


In  this  exploratory  study,  Darcy's  linear  law^  is  used  as  the 
basic  equation  relating  the  appar ent  volumetric  fluid  (air)  flow 
rate  to  the  fluid  pressure  gradient,  i.e.. 


v  u 
a 


i 

r 


U  Sx . 

l 


(2.19) 


where 


a 

i 

u 

r 

i 

v  u 
a  r 


is  air  filled  porosity  of  soil 

is  local  relative  velocity  u1  -  u1 

a  g 

is  the  apparent  flow  velocity  of  air 
is  permeability 
is  the  air  viscosity 


And,  since  Darcy's  law  is  for  steady-state  conditions,  the 
pressure  gradient  on  the  air  is  ]ust  equal  to  the  internal  drag 
torce/volume  on  the  air,  thus 


v  ;)P 

a  a 

9x . 

i 


(2.20) 


‘In  the  general  multipnase 


system 


E  (f 

k-i,K\ 


u 


Drag  formulations  of  the  internal  interaction  force  define 

an  effective  b inale  particle  force,  F1,  which  acts  on  a  single 

4ttR3 

soil  particle  of  radius  R.  Since  — ^ —  is  the  soil  +  air  volume 

s 

per  soil  particle, 


4ttR3 
'  3v 


(2.21) 


In  the  current  study,  a  single  particle  drag  formulation  is 
used  with  three  porosity  regimes.  The  first  regime  applies  to 
relatively  low  porosities,  va^r  less  than  or  equal  to  0.5.  In 
this  regime  the  flow  of  air  through  the  soil  is  modeled  using  a 

g 

drag  formulation  of  Par cy 1 s  law  derived  by  Brinkman  : 


,i  ir,  n  L  R  J.  R  * 

i  =  6,»urR  n  +  7T  +  3k’ 


;  0  <  V  <0.5 

—  3.  — 


(2.22) 


where 


*  viscosity 

=  relative  velocity  of  air  and  soil  particles, 

ui  -  "i 

-  Dp/2  =  characteristic  radius  of  soil  particles 

-  permeability 


In  the  derivation  of  Equation  (2.22)  it  is  assumed  that  the  soil 
particles  are  rigid  spheres  and  that  the  air  flow  is  steady  and 
incompressible.  These  assumptions  may  not  be  adequate  for  a 
final  model  but  they  do  permit  a  first-order  assessment  of  the 
pore-air  phenomena  of  interest. 


The  second  regime  is  transitional  and  applies  for 
0.5  <  v  <  0.8.  In  this  regime  the  soil  is  distended  and  the 
individual  soil  particles  are  not  touching  but  are  still 
relatively  close  together.  The  interaction  between  the  soil  and 
air  is  a  linear  interpolation  on  v  between  Equation  (2.22)  and 
the  single  particle  aerodynamic  drag  interaction  force  [Equation 
(2.24)  below]: 


F 


1 

2 


(0.8  -  v  ) 
a 

0.3 


(va  -  0.5) 

a 

073 


(2.23) 


0.5  l  v  <  0.8 

a 

where  F  is  given  by  Equation  (2.22)  above  and  F3  is  given  by 
Equation  (2.24). 

In  the  third  regime  v  >  0.8,  the  soil  is  highly  distended, 
soil  particles  are  relatively  far  apart  and  do  not  strongly 
interact,  and  the  soil-air  interaction  is  modeled  as  single 
particle  drag: 

Fg  =  !5paffR2CDu^  |u£ 1  ;  va  >  0.8  (2.24) 

where  CD  is  the  variable  drag  coefficient  which,  in  general,  will 
depend  on  the  relative  velocity  and  local  Reynolds  number. 

In  order  to  calculate  the  soil-air  interaction  in  regimes  1 
and  2,  it  is  necessary  to  specify  the  permeability.  The  perme¬ 
ability  of  the  soil  will  change  as  the  soil  is  compressed  by  the 
airblast  and  pores  collapse;  also,  soil  cementation  is  broken  and 
thus  the  particle  size  distribution  is  altered.  It  is  thus 
essential  that  permeability  be  treated  as  a  variable.  As  an 
initial  model,  permeability  is  assumed  to  vary  with  porosity  and 
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characteristic  particle  size.  For  porosities  less  than  or  equal 

to  about  0.48,  an  empirically  based  equation  formulated  by  Carman 
-7 

is  used  . 


R2  45  (1  -  v  ) 2 
a 


;  v  <0.4785 

d  — 


(2.25) 


For  higher  porosities,  a  theoretical  formulation  by  Brinkman  is 
used . 


k  -  1  7  .  4 

R7  "  18  3  +  T - ~  '  3 

a 


-  3)i 


(2.26) 


0.4785  <  v  <  0.8 
a  ~ 


Figure  2.1  shows  the  quantity  k/R  as  a  function  of  porosity  as 
calculated  from  the  Carman  and  Brinkman  equations. 

2.4  TWO-PHASE  FLOW  TEST  CASE 

The  details  of  the  internal  interaction  (or  drag)  forces  and 
the  associated  relative  velocities  (u^)  between  the  air  and  soil 
phases  will  influence  the  wave  propagation  characteristics  in  the 
soil-air  medium.  Under  the  special  case  of  u^  ■  0  which  occurs 
as  the  permeability  goes  to  zero,  the  sound  speed  in  the 
two-phase  medium  is  related  to  the  soil  and  air  densities,  local 
sound  speeds,  and  volume  fractions  as  shown  on  Table  2.1.  Note 
that  when  v  «  v  =0.5,  the  two-phase  sound  speed  reaches  a 

as  g 

minimum  and  is  a  small  fraction  of  the  air  sound  speed  ,  in  fact 


C  "2  ^7^  Ca 


for  v  =  v  1/2 
a  8 


_-.\v 


K/R*»2 


Table  2.1.  Sound  Speed  for  Multiphase  Mixture 


(particles  not  in  contact  with  one  another  and  uj,  = 


1 

'  r2 

p„c  p 

C2  ) 

a  a  / 

v_p,C2  + 

v,peC2 

V„p„  +  v,p. 

sFa  a 

as  s 

sps  aKa 

_ 

- 

.. 

Note  that 


°sC^ 


Thus,  C2 


paCa 

psvsva 


(rigid  soil  particles) 


is  good  approximation  until  v  -*■  0  or  v  -*■  0. 


A  test  case  was  performed  with  the  DICE  code  to  approxi¬ 
mately  simulate  the  specialized  condition  just  described  and  to 
verify  the  ability  to  calculate  low  wave  speeds  in  two-phase 
mixtures.  The  problem  configuration  and  soil/air  volume 
fractions  are  indicated  below: 


1  bar 

overpressure 

shock 

conditions 


Air  Only 


Soil -Air  Mixture 


Distance  (m) 


A  1  bar  overpressure  shock  wave  travels  through  100  m  of  air  at 
STP  in  a  simulated  shock  tube  and  then  interacts  with  the 
soil/air  mixture. 


Figure  2.2  shows  the  overpressure  versus  distance  profiles 
at  three  times  prior  to  the  shock  wave  interaction  with  the 
two-phase  mixture.  Figure  2.3  shows  the  overpressure  profiles  at 
four  times  after  the  shock  wave  begins  interacting  with  the 
soil/air.  Figure  2.4  is  the  related  shock  speed  versus  time 
curve.  The  shock  wave  in  the  air  propagates  at  ~463  m/s;  when 
the  initial  shock  reaches  the  relative  dense  soil/air  mixture 
(p  =  1  gm/cm3)  at  x  *  100  m,  a  reflected  shock  wave  travels  back 
into  the  air  with  a  reflection  factor  of  ~2.6. 
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Figure  2.2.  Overpressure  Versus  Distance  Profiles  at  t  =  .096,  .166  and  .208  Seconds 
for  Shock  Tube  Test  Case. 


■  r  r 


Figure  2.3.  Overpressure  Versus  Distance  Profiles  at  t 
.448  Seconds  for  Shock  Tube  Test  Case. 
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2.5  AIR  ONLY  TEST  CASE  FOR  WEAK  SHOCK  (AP  -  0.1  bar) 


For  some  airblast  problems  of  interest,  the  air  overpressures 
are  quite  low.  Therefore,  a  low  overpressure  test  integration  of 
the  DICE  code  was  performed  to  demonstrate  that  the  code  can 
calculate  shocks  and  rarefaction  waves  correctly  at  overpressures 
of  ~0.1  bar.  The  following  sketch  shows  the  initial  pressure 
conditions  used: 


Initially,  the  air  was  at  rest  in  both  regions  and  is  assumed 
to  follow  the  ideal  gas  law  (7  *  1.4).  The  zoning  consisted  of  a 
uniform  grid  with  one  meter  cells. 

Analytical  (closed  form)  solutions  (see  Ref.  10,  pp.  181-191) 
can  be  obtained  in  1-D  by  considering  the  equations  governing  the 
conservation  of  mass,  momentum  and  energy  plus  the  ideal  gas 
equation  of  state  for  air.  The  results  for  the  initial  conditions 
considered  are  summarized  in  Figure  2.5.  In  this  x-t  diagram,  the 
theoretical  positions  of  the  shock  front,  contact  discontinuity 
and  rarefaction  fan  are  plotted  as  a  function  of  time.  The 
physical  parameters  describing  the  properties  of  each  of  the 
separate  regions  are  also  shown. 


Figure  2.5.  Theoretical  Positions  of  Shock  Front,  Contact  Discontinuity  and 
Rarefaction  Fan  Versus  Time  for  Low  Overpressure  Test  Case. 


DICE  code  integration  results  at  .030  s,  as  well  as  the 
theoretical  solution,  are  shown  in  Figures  2.6  through  2.9.  Air 
overpressure  versus  distance  is  shown  in  Figure  2.6.  A  shock  wave 
of  0.1  bar  overpressure  is  seen  to  move  forward  into  the  undis¬ 
turbed  region,  while  rarefaction  waves  move  backward  into  the 
initial  overpressure  region.  These  rarefaction  waves  comprise  the 
rarefaction  fan.  For  both  the  shock  front  and  rarefaction  fan, 
DICE  is  able  to  define  the  rapid  variations  with  five  cells.  The 
shock  pressure  is  matched  almost  exactly.  Similar  accurate 
results  are  seen  for  the  discontinuities  within  the  velocity, 
temperature  and  density  profiles  shown  in  Figures  2.7,  2.8  and 
2.9,  respectively.  We  conclude  that  DICE  can  handle  shocks  and 
rarefaction  waves  correctly  for  low  overpressures  at  the  0.1  bar 
level . 
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Figure  2.6.  Air  Overpressure  Versus  Position  at  t  =  .030  Seconds  for  Low 
Overpressure  Test  Case. 
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Figure  2.8.  Air  Temperature  Versus  Position  at  t  =  .030  Seconds  for  Low 
Overpressure  Test  Case. 
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gure  2.9.  Air  Density  Versus  Position  at  t  =  .030  Seconds  for  Low  Overpressure 
Test  Case. 


SECTION  3 


PORE -AIR  PHENOMENA  IN  A  SHOCK  TUBE 

3.1  LABORATORY  EXPERIMENTS 


3 

Pore-air  phenomena  were  studied  by  Zernow  et  al.  in  a 
series  of  laboratory  experiments.  The  experiments  involved  the 
use  of  a  shock  tube  (Figure  3.1)  to  load  the  surface  of  a  column 
of  sand  with  a  known  overpressure  history.  By  breaking 
diaphragms  1  and  2  at  appropriate  times,  the  surface  of  the 
column  of  sand  was  exposed  to  air  overpressure  histories  similar 
to  the  positive  phase  of  the  airblast  at  the  100  psi  range  from  a 
1  Mt  nuclear  burst.  No  negative  phase  was  simulated.  The  motion 
of  the  sand  column  was  observed  photographically  and  pore-air 
pressure  measurements  were  made  at  various  depths  in  the  sand. 

Figure  3.2  3hows  the  observed  surface  load  overpressure  time 
histories  for  the  experiments  conducted  with  a  8"  diameter  shock 
tube.  Two  orifice  diameters  were  used  to  control  the  overpres¬ 
sure  decay  rate.  The  1.78"  diameter  orifice  yielded  a  faster 
decay  rate  with  a  shorter  duration  (-.450  s  positive  phase) 
compared  to  1.36"  diameter  orifice  with  a  positive  phase  of 
-.750  s.  Also  shown  in  the  figure  is  an  overpressure  time 
history  described  by  a  simple  exponential  decay  function  of  the 
form 
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Figure  3.2.  Overpressure  Load  to  Surface  of  Sand  Column  Versus  Time 
Observed  in  Experiments  (Zernow  et  al .,  1973)  and  Fit 
Exponentially  (used  in  DICE  Cases  1  to  4  Calculations). 


where 


APT  is  the  overpressure  at  time  t 

Li 

APq  is  the  peak  overpressure  (at  the  shock  front) 

t  is  the  elapsed  time  (after  peak  overpressure  arrival) 
x  is  the  decay  parameter  constant. 

A  value  of  6.89  bars  (100  psi)  for  APq  and  .174  s  for  x 
results  in  an  approximate  fit  to  both  of  the  experimentally 
observed  time  histories.  This  loading  function  is  used  in 
subsequent  DICE  code  numerical  simulations. 

The  average  density  of  the  standard  sand  used  was  found  to 
be  1.60  gm/cm3 .  The  estimated  grain  density  of  2.40  gm/cm3 
corresponds  to  a  porosity  (or  air  void  ratio)  of  0.33.  Figure 
3.3  shows  the  particle  size  distribution  (%  by  weight)  for  the 
standard  sand;  the  mass  mean  diameter  is  -650  n.  Some  of  the 
experiments  were  conducted  with  "coarse"  sand  with  particle 
diameters  greater  than  841  n  sieved  from  the  standard  sand.  The 
coarse  sand  had  an  average  density  of  1.55  gm/cm 3  and  porosity  of 
0.354. 


The  bulk  compressibility  modulus  was  determined  to  be 
85,000  psi  (-5860  bars),  which  implies  a  compression  wave 
velocity  of  -605  m/s.  For  the  6  ft  sand  column  this  means  that 
the  transit  time  for  the  initial  shock  to  reach  the  bottom  is 
.003  s. 

Note  that  in  the  experiments,  the  bottom  of  the  sand  column 
(-6  ft)  was  restricted  from  moving  downward.  Ideally,  the  sand 
column  should  be  "infinitely  long"  with  respect  to  the  physical 
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effects  which  are  being  observed.  In  an  ideal  elastic  material, 
the  sand  at  a  depth  equivalent  to  the  bottom  of  the  sand  column 
would  behave  similar  to  the  sand  near  the  surface.  A  pressure 
force  applied  to  the  sand  surface  would  be  felt  3  ms  later  at  a 
6  ft  depth,  and  would  eventually  cause  the  entire  sand  column  to 
translate  downward.  However,  the  rigid  bottom  in  the  experiments 
prevents  this  downward  translation  at  the  surface  as  well  as  at 
the  bottom  since  the  sand  is  essentially  elastic.  To  simulate  an 
infinitely  long  sand  column  with  respect  to  pore-air  phemonena, 
the  bottom  closure  of  the  shock  tube  had  small  holes  in  it  to 
allow  restricted  air  flow.  The  intent  was  to  approximate  a 
permeability  similar  to  that  of  the  sand  above  the  bottom 
closure. 

Figure  3.4  shows  the  vertical  displacement  of  the  sand 
surface  as  a  function  of  time  from  the  1.36"  and  1.78"  diameter 
orifice  tests.  The  initial  small  displacements  are  considered  to 
be  a  "precursor"  sand  cloud  resulting  from  the  shock  tube  induced 
"initial  large  amplitude  reverberation  oscillations  in  the  very 
early  portion  of  the  shock  pulse"  affecting  the  top  -2  cm  of  the 
sand  column.  The  main  sand  column  starts  to  move  upward  at 
t  -.400  s  for  the  standard  sand  1.36"  tests  and  maintains  a  nearly 
constant  velocity  of  ~43  cm/s  for  another  .300  to  .400  s.  The 
1.78"  tests  with  faster  overpressure  decay  rates  resulted  in  the 
main  sand  column  moving  upward  at  an  earlier  time  (-.200  s)  and 
higher  velocity  (~77  cm/s).  Both  sets  of  curves  for  the  standard 
sand  experiments  indicate  that  the  maximum  height  attained  by  the 
lofted  sand  is  nearly  20  cm.  If  the  only  force  acting  on  a  sand 
particle  projected  upward  were  gravity,  the  initial  lofting 
velocity  needed  to  reach  20  cm  height  is  V  *  /2gH  -200  cm/s. 
Including  air  drag  effects  would  increase  the  needed  initial 
lofting  velocity.  The  observed  maximum  velocity  of  80  cm/s  would 
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Figure  3.4.  Main  Surface  Displacement  Versus  Time  for  8"  Diameter 
Shock  Tube  Experiments  (Zernow  et  at.,  1973). 


project  a  sand  particle  less  than  4  cm  above  the  surface.  Thus, 
the  experimental  evidence  points  to  the  sustained  application  of 
an  additional  upward  force  associated  with  pore-air  phenomena. 

Figure  3.4  also  shows  the  main  surface  displacement  for 
experiments  run  with  the  coarse  sand  (particle  diameter  >841  m) . 
The  coarser  sand  particles  yielded  an  increase  in  porosity  and 
permeability.  And  although  the  same  1.36"  diameter  orifice 
pressure  loading  was  used,  there  was  a  reduction  in  the  main 
column  motion  by  a  factor  of  4  both  in  velocity  and  displacement 
when  compared  to  the  standard  sand  cases.  The  coarse  sand 
characteristics  are  very  uncertain;  a  possible  explanation  is 
described  in  Section  3.3. 

Figure  3.5  shows  the  position  of  the  near-surface  marked 
sand  layers  as  a  function  of  time  for  the  1.78"  and  1.36" 
diameter  orifice  tests.  Sustained  upward  motion  is  again  seen 
for  the  sand  below  the  surface  to  depths  of  approximately 
10  to  15  cm. 

Additional  experiments  involved  saturating  the  standard  sand 
100%  with  water  so  that  the  porosity  and  permeability  were 
essentially  zero,  and  saturating  the  standard  sand  to  50% 
saturation  level.  The  result  was  either  total  or  near  total 
elimination  of  main  column  motion. 

3.2  NUMERICAL  SIMULATIONS 

Numerical  simulations  were  performed  with  the  multiphase 
DICE  code  in  order  to  examine  the  pore-air  phenomenology  in  a 
homogeneous  porous  soil.  The  calculations  were  designed  using 
one-dimensional  (1-D)  planar  geometry  to  approximate  the 
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laboratory  shock  tube  experiments.  The  equations  governing  the 
relevant  physics  are  discussed  in  Section  2. 

Figure  3.6  defines  the  physical  parameters  used  in  DICE 
Cases  1  to  4 .  The  overpressure  loading  function  applied  at  the 
soil  surface  is  the  exponential  fit  shown  in  Figure  3.2,  falling 
between  the  two  overpressure  time  histories  measured  in  the 
laboratory  experiments.  Its  rate  of  decay  is  slower  than  for  the 
1.78"  diameter  orifice  tests  but  faster  than  for  the  1.36"  tests. 


The  soil  model  parameters  are  taken  from  Reference  3.  The 
initial  soil  density  is  1.6  gm/cm3  with  a  porosity  of  v  =  0.33. 
The  soil  is  treated  as  a  hydro-elastic  material.  The  605  m/s 
compression  wave  velocity  implies  a  3.3  ms  transit  time  for  the 
initial  shock  pulse  at  the  surface  to  reach  the  bottom  grid 
boundary.  A  semi-inf initely  long  column  is  modeled  by  imposing  a 
continuum  boundary  at  the  grid  bottom,  allowing  soil  and  air  to 
flow  past  the  bottom;  note  that  in  the  laboratory  experiments, 
the  soil  is  rigidly  constrained  at  the  bottom  of  the  shock  tube. 


In  the  calculation  the  soil  is  modeled  using  only  one 
characteristic  particle  3ize.  Case  1  uses  a  characteristic 
particle  size  diameter  of  =  0.1  cm  (all  particles  assumed  to 
be  uniform  spheres),  while  Case  2  uses  =  0.05  cm.  The 
permeability  model  described  in  Section  2.3  indicates  that  the 
permeability  increases  with  porosity,  and  for  a  given  porosity, 
the  permeability  varies  as  the  square  of  the  particle  size. 

Thus,  the  permeabiity  of  the  soil  in  Case  1  is  four  times  larger 
than  in  Case  2  for  a  given  porosity.  Cases  1  and  2  include 
effects  of  a  uniform  gravitational  field  (g  =  980  cm/s2). 
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Figure  3.6.  Pressure  Loading  and  Soil  Parameters  for  1-D  Shock 
Tube  Calculations  (Cases  1  to  4). 


Cases  2  and  3  compare  the  effect  of  gravity  on  the  pore-air 
lofting  process.  =  0.05  cm  for  each  case;  but  in  Case  3,  zero 

gravity  is  assumed.  Case  4,  which  also  has  zero  gravity,  has  a 
factor  of  ~3  increase  in  particle  size  and  a  factor  of  ~8 
increase  in  permeability  from  Case  3. 

Figure  3.7  shows  the  1-D  grid  used  in  the  calculations.  The 
surface  soil  cell  is  2  cm,  with  each  subsequent  cell  10%  larger 
than  the  cell  above  it.  All  cells  above  the  original  surface 
(z  -  0)  are  of  uniform  size  (2  cm). 

3.2.1  Case  1  (D^  =  0.1  cm,  g  -  980  cm/s2) 

Figure  3.8  shows  the  initial  pressure  distribution  versus 
depth  for  soil  and  air  for  Case  1.  The  total  stress  shown  is  the 
sum  of  the  pressure  due  to  the  pore  air  and  the  effective  stress 
contribution  from  the  soil.  The  initial  air  pressure  was  taken 
to  be  1.013  bars  (=  1  atm)  throughout  the  soil  column  depth.  The 
initial  lithostatic  pressure  at  depth  in  the  soil  column  is  the 
pressure  required  to  support  the  weight  of  the  total  amount  of 
soil  above  that  depth.  It  increases  essentially  linearly  from 
zero  at  the  surface  to  a  value  of  0.314  bars  at  a  depth  of 
200  cm. 


Figure  3.9  shows  the  pressure  and  velocity  distributions  for 
soil  and  air  at  t  -  .001  s  after  the  initial  application  of  the 
pressure  loading  function  at  the  soil  surface.  The  soil 
compression  wave  velocity  of  605  m/s  implies  a  transit  distance 
of  60  cm  by  this  time  for  the  initial  pressure  loading  impulse. 
First  principle  shock  theory  for  an  elastic  medium  yields  a  soil 
particle  velocity  of  71  cm/s  behind  the  shock  front  for  an 
initial  pressure  of  6.89  bars.  DICE  gives  good  agreement  with 
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DEPTH  (cm)  1. 008  MSEC 


Figure  3.9.  (a)  Soil,  Air  and  Total  Stress  Versus  Depth,  and 

(b)  Soil  and  Air  Velocity  Versus  Depth  at 
t  =  .001  Seconds  for  Case  1  (D  =0.1  cm). 


these  theoretical  values.  The  pressure  loading  shock  front  is 
seen  to  travel  to  the  theoretical  depth  (the  numerically  smeared 
transition  region  across  the  shock  front  is  centered  at  ~6Q  cm); 
the  downward  soil  velocity  behind  the  front  is  approximately 
70  cm/s. 

At  t  -  .001  s,  the  soil  is  supporting  all  of  the  overpres¬ 
sure  load  except  near  the  surface,  where  air  traveling  down  at 
over  1500  cm/s  has  begun  to  permeate  into  the  soil.  As  the  air 
is  forced  down  into  the  soil  matrix,  the  pore  air  density  and 
thus  pore  air  pressure  increases  and  the  air  supports  a  portion 
of  the  pressure  load. 

Figure  3.10  shows  that  by  t  »  .050  s  the  pressure  loading  of 
the  soil  surface  has  affected  the  entire  column.  Below  100  cm 
the  entire  overpressure  load  is  supported  by  the  soil  (the  air  is 
still  at  its  initial  value).  Above  this  depth  the  pore  air 
pressure  increases  nearly  linearly  to  the  surface.  At  the 
surface  the  pressure  load  is  supported  almost  entirely  by  the 
air.  At  the  surface  the  air  is  flowing  down  through  the  soil  at 
-200  cm/s.  The  peak  downward  air  velocity  of  nearly  350  cm/s  is 
at  a  depth  of  ~50  cm.  Below  100  cm  the  soil  and  air  matrix  move 
downward  at  a  velocity  54  cm/s  consistent  with  the  theoretical 
velocity  behind  a  shock  front  of  overpressure  AP^(t  “  .050)  • 
6.89e_L/-i74  =  5.17  bars. 

The  total  pressure  at  any  time  will  still  contain  the 
initial  air  and  soil  contributions  discussed  previously  (see 
Figure  3.8).  In  addition,  because  of  the  decay  in  time  of  the 
overpressure  loading  function,  there  is  a  net  overpressure 
differential  from  the  surface  to  any  depth,  due  to  the  transit 
time  (at  sound  speed  velocity  C)  required  for  information  (the 
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Figure  3.10.  (a)  Soil,  Air  and  Total  Stress  Versus  Depth,  and 

(b)  Soil  and  Air  Velocity  Versus  Depth  at 
t  =  .050  Seconds  for  Case  1  (D  =0.1  cm). 


pressure  load)  from  the  surface  to  get  to  that  depth.  The 
theoretical  overpressure  at  depth  z  and  time  t  due  to  the  loading 
function  is  thus  AP(z,t)  =  6.89e  ^  z/c)/*174  ^ars.  The 
overpressure  differential  Ap  from  the  surface  to  depth  z  (as  a 
function  of  time)  is  thus 

Ap(z,t)  -  6.89e-t/*174[eZ/C/*174  -1]  bars.  (3.2) 

It  is  this  net  overpressure  differential  which  provides  the 
upward  force  which  accelerates  the  soil  and/or  air  upward,  thus 
decreasing  its  downward  velocity.  Note  the  exponential  decay  in 
time  which  will  result  in  weaker  upward  forces. 

Figure  3.11  shows  that  by  t  -  .100  s  the  air  has  begun  to 
flow  through  the  soil  at  a  depth  of  over  150  cm,  with  the  peak 
downward  velocity  of  ~200  cm/s  at  ~80  cm  depth.  More 
significantly,  by  t  -  .100  s  Ulfe  Air  is  supporting  the  entire 
presure  load  at.  the  surface.  The  soil  particles  become 
fluidized;  they  are  no  longer  able  to  support  any  loads.  The 
soil  particle  motion  in  a  fluidized  state  is  dominated  by  gravity 
and  internal  drag/permeability  interactions  with  the  air. 

Figure  3.12  shows  that  at  t  =  .200  s  the  soil  is  fluidized 
to  a  depth  of  ~40  cm  and  a  pressure  gradient  causing  upward 
accelerations  exists  in  this  near-surface  region.  Much  lower 
sound  speeds  in  the  fluidized  soil  as  compared  to  the  unfluidized 
soil  (see  Section  2.4)  will  result  in  larger  upward  pressure 
forces  as  indicated  by  Equation  (3.2).  And  because  the  air 
density  is  about  three  orders  of  magnitude  less  than  the  soil 
density,  there  is  a  larger  increase  in  upward  acceleration  of  the 
air  as  compared  to  the  soil  particles.  The  results  of  this  are 
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Figure  3.12.  (a)  Soil,  Air  and  Total  Stress  Versus  Depth,  and 

(b)  Soil  and  Air  Velocity  Versus  Depth  at 
t  =  .200  Seconds  for  Case  1  (D  =0.1  cm). 


seen  in  the  velocity  profiles  in  Figure  3.12.  The  air  in  the 
fluidized  region  is  moving  upward  relative  to  the  soil.  Internal 
force  interactions  with  the  soil  particles  is  providing  a 
deceleration  of  the  air,  while  causing  the  soil  particles  to 
accelerate  upward;  the  soil  in  the  fluidized  region  is  in  fact 
moving  upward  relative  to  the  remainder  of  the  soil  column. 

Figure  3.13  shows  that  by  t  -  .400  s  the  region  of 
fluidization  extends  to  ~115  cm  depth.  The  soil  column  above  100 
cm  depth  is  moving  upward,  with  a  maximum  velocity  of  67  cm/s  at 
the  surface.  The  upward  air  velocity  is  also  largest  at  the 
surface,  with  a  value  of  over  100  cm/s. 

Figure  3.14  shows  that  at  t  =  .600  s  the  region  of 
fluidization  extends  below  150  cm  depth.  The  soil  column  above 
this  depth  is  moving  upward,  with  a  maximum  velocity  of 
40  cm/s  near  the  surface.  The  upward  air  velocity  peaks  at  the 
surface  with  a  value  of  over  80  cm/s . 

The  progression  of  fluidization  down  the  soil  column  can  be 
seen  in  the  soil  density  depth  profiles  at  various  times  on 
Figure  3.15.  Up  to  t  »  .100  s  the  soil  density  is  still 
essentially  at  its  initial  value  of  1.6  gm/cm3 .  The  region  of 
fluidization  is  clearly  seen  in  the  profiles  at  later  times.  The 
soil  density  drops  below  1.6  gm/cm3  as  the  soil  becomes 
increasingly  distended.  By  t  »  .600  s  porosities  near  the 
surface  have  increased  to  over  0.5. 

During  the  numerical  simulations  various  predicted  physical 
quantities  are  monitored  at  several  station  depths  in  the  1-D 
grid.  This  allows  for  time  histories  of  quantities  such  as  air 
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Figure  3.13.  (a)  Soil,  Air  and  Total  Stress  Versus  Depth,  and 

(b)  Soil  and  Air  Velocity  Versus  Depth  at 
t  =  .400  Seconds  for  Case  1  (D  =0.1  cm). 
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Figure  3.15.  Soil  Density  Versus  Depth  at  Varioi 


overpressure  and  soil  velocity  at  the  selected  station  depths  to 
be  analyzed. 

Figure  3.16  shows  the  air  overpressure  versus  time  for  Case 
1  at  various  station  depths  of  the  modeled  soil  column  from  5  cm 
to  175  cm.  The  peak  overpressure  decreases  with  depth  while  the 
time  of  arrival  of  the  peak  overpressure  increases  with  depth. 
This  figure  shows  that  by  t  =  .200  s,  air  overpressures  down  to 
40  cm  depth  have  become  approximately  equivalent  to  the  over¬ 
pressure  loading  applied  at  the  surface.  Thus,  above  40  cm  depth 
the  soil  has  become  fluidized  and  is  supported  entirely  by  the 
air  . 


Figure  3.17  shows  the  soil  velocity  versus  time  at  various 
station  depths  from  5  cm  to  175  cm.  Up  to  time  t  -.100  s  the 
soil  column  essentially  moves  downward  as  a  rigid  body. 
Fluidization  at  each  station  depth  becomes  apparent  after  this 
time  when  the  soil  flowing  at  that  depth  increases  its  velocity 
above  the  velocity  of  the  soil  column  below  that  depth.  At  the 
5  cm  depth  station  a  maximum  upward  velocity  of  -65  cm/s  is 
achieved  at  t  -.400  s,  after  which  time  gravity  becomes  the 
dominating  force. 

3.2.2  Case  2  (D  0.05  cm,  g  =*  980  cm/s2) 

P 

The  effect  of  permeability  of  the  modeled  soil  is  examined 
in  Case  2  by  using  a  characteristic  particle  size  diameter  of 
D  =  0.05  cm,  which  is  half  that  used  in  Case  1.  For  a  qiven 
porosity,  the  soil  permeability  in  Case  2  is  four  times  smaller 
th._..n  in  Case  1  (see  the  permeability  model  described  in  Section 
2.3). 
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Air  Overpressure  Versus  Time  at  Various  Station 
Depths  for  Case  1  (D  =0.1  cm). 
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Figure  3.18  shows  the  pressure  and  velocity  distributions 
for  soil  and  air  at  t  *  .001  s  for  Case  2.  Air  is  flowing  down 
through  the  soil  at  just  over  200  cm/s  for  Case  2  as  compared  to 
over  1500  cm/s  for  Case  1.  There  is  also  a  corresponding 
decrease  in  the  amount  of  support  of  the  overpressure  load  at  the 
surface  by  the  air. 

By  t  =  .050  s  (Figure  3.19)  the  peak  downward  velocity  of 
-200  cm/s  is  found  at  a  depth  of  25  cm  for  Case  2,  slower  and 
shallower  than  the  -330  cm/s  peak  at  50  cm  depth  for  Case  1.  The 
air  pressure  also  indicates  that  there  is  a  decrease  in  permea¬ 
tion  of  the  soil;  air  is  supporting  at  least  part  of  the  pressure 
loading  down  to  almost  100  cm  depth  for  Case  1,  but  only  to  50  cm 
for  Case  2. 

Fluidization  at  the  surface  occurs  near  t  =  .100  s  (Figure 
3.20)  for  both  Case  1  and  Case  2,  when  the  air  is  supporting  the 
entire  pressure  load  at  the  surface.  By  t  *  .200  s,  Figure  3.21 
shows  the  soil  fluidized  to  only  ~25  cm  depth  for  Case  2,  com¬ 
pared  to  -40  cm  depth  for  Case  1  (Figure  3.12).  By  t  =  .400  s 
(Figure  3.22)  the  soil  in  Case  2  is  fluidized  to  a  depth  of 
~70  cm  compared  to  ~120  cm  for  Case  1,  indicating  the  effect  of 
the  decrease  in  permeability  for  Case  2. 

The  air  also  achieves  smaller  upward  velocities  compared 
with  those  of  Case  1.  This  is  seen  in  Figure  3.22  at  t  =  .400  s 
and  Figure  3.23  at  t  =  .600  s.  Evidence  of  the  increased 
air/soil  drag  forces  due  to  the  decrease  in  permeability  is  found 
in  the  smaller  relative  velocities  between  the  air  and  soil 
particles  in  the  fluidized  region. 
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Figure  3.19. 


25  50  75  100  125  150  175  200 

DEPTH  (cm)  50.  00  MSEC 

(a)  Soil,  Air  and  Total  Stress  Versus  Depth,  and 

(b)  Soil  and  Air  Velocity  Versus  Depth  at 

t  =  .050  Seconds  for  Case  2  (D  =0.05  cm). 
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(b)  Soil  and  Air  Velocity  Versus  Depth  at 
t  =  .200  Seconds  for  Case  2  (D  =0.05  cm). 
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Figure  3.23.  (a)  Soil,  Air  and  Total  Stress  Versus  Depth,  and 

(b)  Soil  and  Air  Velocity  Versus  Depth  at 

t  =  .600  Seconds  for  Case  2  (D  =  0.05  cm). 


Figure  3.24  shows  the  air  overpressure  time  histories  at 
various  station  depths  far  Case  2.  The  peak  overpressure  at  a 
particular  depth  is  smaller  than  in  Case  1  and  also  arrives 
later . 

Figure  3.25  shows  the  soil  velocity  time  histories  at 
various  station  depths  for  Case  2.  Fluidization  is  indicated  by 
the  increase  in  velocity  at  a  particular  depth  above  the  velocity 
of  the  soil  column  below  that  depth.  Fluidization  at  a  partic¬ 
ular  depth  below  the  surface  occurs  at  later  times  for  Case  2 
than  for  Case  1.  Equation  3.2  indicates  that  the  later  time 
leads  to  a  weaker  upward  pressure  force  at  that  depth  for  the 
fluidized  soil  in  Case  2.  At  the  5  cm  depth  station  a  maximum 
upward  velocity  of  ~35  cm/s  is  achieved  at  t  -.400  s.  This  peak 
velocity  is  less  than  in  Case  1  by  about  a  factor  of  two. 

However,  the  time  at  which  the  peak  occurs  is  approximately  the 
same  in  both  cases . 

3.2.3  Case  3  (D  =  0.05  cm,  g  =  0) 

P 

The  effect  of  gravity  on  the  modeled  soil  column  is  examined 

in  Case  3  by  setting  g  =  0  instead  of  980  cm/s2  as  used  in  Cases 

1  and  2.  A  characteristic  particle  size  of  diameter  D  =  0.05  cm 

IP 

is  used  so  that  this  case  compares  directly  to  Case  2.  In  the 
absence  of  gravity  there  is  no  initial  lithostatic  pressure  at 
depth  in  the  soil  column  required  to  support  the  soil  weight 
above  that  depth.  There  is  only  an  initial  pore  air  pressure  of 
1.013  bars  (-  1  atm)  taken  throughout  the  soil  column  depth. 

Figures  3.26  through  3.31  show  the  pressure  and  velocity 
distributions  for  soil  and  air  from  times  t  =  .001  s  to 
t  =  .600  s  for  Case  3.  Op  to  t  -  .100  s  the  air  pressure  and  the 
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Figure  3.25.  Soil  Velocity  Versus  Time  At  Various  Station  Depths 
for  Case  2  (D  =  0.05  cm). 
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Figure  3.26.  (a)  Soil,  Air  and  Total  Stress  Versus  Depth,  and 

(b)  Soil  and  Air  Velocity  Versus  Depth  at 
t  =  .001  Seconds  for  Case  3  (D  =0.05  cm). 
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Figure  3.27.  (a)  Soil,  Air  and  Total  Stress  Versus  Depth,  and 

(b)  Soil  and  Air  Velocity  Versus  Depth  at 
t  =  .050  Seconds  for  Case  3  (D  =0.05  cm). 
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(a)  Soil,  Air  and  Total  Stress  Versus  Depth,  and 

(b)  Soil  and  Air  Velocity  Versus  Depth  at 

t  =  .200  Seconds  for  Case  3  (D  =0.05  cm). 
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soil  and  air  velocity  distributions  are  essentially  the  same  as 
those  of  Case  2.  The  soil  pressure  differences  reflect  the 
initial  lithostatic  pressure  field  imposed  in  Case  2. 

By  t  -  .200  s  the  soil  is  fluidized  to  ~25  cm  depth  for  both 
Case  3  and  Case  2.  For  Case  3  the  relative  velocity  between  the 
air  and  soil  in  this  region  is  almost  zero,  the  air  being 
slightly  larger  in  the  upward  direction.  The  relative  velocity 
is  larger  for  Case  2  due  to  gravity,  which  accelerates  the  soil 
downward  in  the  opposite  direction  from  the  upward  accelerating 
air.  Internal  interaction  forces  will  influence  the  motion  until 
there  is  no  relative  velocity  between  the  soil  and  air. 

By  t  -  .400  s,  the  soil  is  fluidized  to  ~75  cm  depth  for 
Case  3  compared  to  --70  cm  depth  for  Case  2.  The  difference 
becomes  larger  with  increasing  depth  because  for  Case  2,  the  air 
pressure  at  a  particular  depth  must  overcome  the  initial  litho¬ 
static  pressure  as  well  as  the  pressure  in  the  soil  resulting 
from  the  surface  pressure  loading  before  the  soil  becomes 
fluidized. 

The  upward  velocity  of  the  air  in  the  fluidized  region  is 
approximately  the  same  for  both  cases  at  t  -  .400  s  and 
t  -  .600  s.  But  again,  the  relative  velocity  is  larger  for  Case 
2  as  gravity  acts  on  the  soil.  By  t  -  .600  s  the  soil  particles 
are  in  fact  moving  upward  faster  than  the  air  for  Case  3. 

Figure  3.32  shows  the  air  overpressure  time  histories  at 
various  station  depths  for  Case  3.  When  compared  to  Case  2 
j  (Figure  3.23)  they  appear  almost  identical.  However,  for  Case  3 

the  overpressures  all  fall  along  the  same  decay  path  once  the 
soil  has  been  fluidized  at  the  various  depths.  In  Case  2,  each 
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successively  deeper  station  required  a  larger  air  overpressure  to 
fluidize  the  soil,  due  to  the  initial  lithostatic  field. 

Figure  3.33  shows  the  soil  velocity  time  histories  at 
various  station  depths  for  Case  3.  The  fluidized  soil  velocities 
are  larger  than  in  Case  2,  again  due  to  the  lack  of  gravity  in 
Case  3.  The  peak  soil  velocity  for  the  5  cm  station  depth  is 
-45  cm/s  in  Case  3  compared  to  -35  cm/s  in  Case  2. 

3.2.4  Case  4  (D^  =  0.142  cm,  g  -  0) 

The  effect  of  increased  permeability  of  the  soil  column 

without  the  effects  of  gravity  is  examined  in  Case  4  by  using  a 

characteristic  particle  size  diameter  of  D  -  0.142  cm.  For  a 

P 

given  porosity,  the  permeability  model  described  in  Section  2 
yields  a  soil  permeability  in  Case  4  which  is  over  eight  times 
larger  than  in  Case  3. 

Figures  3.34  through  3.39  show  the  pressure  and  velocity 
distributions  for  soil  and  air  from  time  t  *  .001  s  to  t  *  .400  s 
for  Case  4.  The  increase  in  permeability  is  seen  immediately  at 
t  »  .001  sec.  Air  is  flowing  down  through  the  soil  at  almost 
5000  cm/s  compared  to  -200  cm/s  for  Case  3.  Air  is  supporting  at 
least  part  of  the  pressure  loading  down  to  -20  cm  depth  for  Case 
4,  compared  to  -5  cm  for  Case  3. 

By  t  =  .050  s  (Figure  3.35)  the  peak  downward  velocity  of 
-450  cm/s  is  at  a  depth  of  75  cm  for  Case  4,  faster  and  deeper 
than  the  -200  cm/s  peak  at  25  cm  depth  for  Case  3.  The  air 
pressure  also  indicates  the  increase  in  permeation  of  the  soil; 
the  air  is  supporting  at  least  part  of  the  pressure  loading  down 
to  over  125  cm  depth  for  Case  4,  but  only  to  50  cm  for  Case  3. 


VEL3CITY  (CM/S) 


VELOCITY  (CM/SEC) 


CALIFORNIA  RESEARCH  AND  TECHNOLOGY.  INC. 
x  .  DICE  3280  CASE  4  (DP=. 142  CM.  G-OJ  PRESSURES  VS.  DEPTH 

)  10*  j  "1  T  I  ~f - ]  - 1  •  I  ' 


in 

cc 

ct 

cn  o. 


cc 

=>  . 
<n  4. 
tn 

LU 

cc 

a. 


Initial  Pressure  Loading 


■Total 


0  25  50 


75  100 


’  *  1 2  5  *  *  *  *1 50  "  *1*7  5  ^  *  *2 


DICE  3280  CASE  4  (QP-. 142  CM.  G-OJ  VELOCITIES  VS. . DtPTH 


-1000 


-2000 


-3000 


i«— Ail 


:  i 

-4000  t-  j 

....!  / 


-5000 


-6000  1 - - — 1 

-25  0 


25  50 


75 

DEPTH  (cm) 


1. 008  MSEC 


Figure  3.34.  (a)  Soil,  Air  and  Total  Stress  Versus  Depth,  and 

(b)  Soil  and  Air  Velocity  Versus  Depth  at 

t  =  .001  Seconds  for  Case  4  (D  =0.142  cm). 
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Figure  3.35.  (a)  Soil,  Air  and  Total  Stress  Versus  Depth,  and 

(b)  Soil  and  Air  Velocity  Versus  Depth  at 

t  =  .050  Seconds  for  Case  4  (D  =0.142  cm). 
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(b)  Soil  and  Air  Velocity  Versus  Depth  at 

t  =  .200  Seconds  for  Case  4  (D  =0.142  cm). 


VELOCITY  (CM/SEC  1 


DEPTH  (cm)  0.400  SEC  >Vv 

V«Vv 

Figure  3.39.  (a)  Soil,  Air  and  Total  Stress  Versus  Depth,  and 

(b)  Soil  and  Air  Velocity  Versus  Depth  at 

t  =  .400  Seconds  for  Case  4  (D  =0.142  cm). 

r  •.*••** 


> 


107 


By  t  -  .100  8  pore  air  pressure  has  increased  to  over  175  cm 
depth  compared  to  ~75  cm  depth  for  Case  3. 

Figure  3.37  shows  that  by  t  =  .200  s  fluidization  of  the 
soil  has  occurred  to  ~60  cm  depth  for  Case  4  compared  to  only 
~2S  cm  depth  for  Case  3.  By  t  -  .300  s  Figure  3.38  shows 
fluidization  down  to  almost  140  cm  depth. 

At  t  =  .400  s  the  entire  200  cm  depth  column  is  fluidized  in 
Case  4,  compared  to  ~70  cm  for  Case  3.  Figure  3.39  shows  that 
air  is  supporting  the  entire  load;  the  entire  fluidized  column  is 
moving  upward,  with  peak  velocity  at  the  surface  of  -190  cm/s  for 
the  air  and  almost  150  cm/s  for  the  soil  particles.  For  Case  3 
both  the  air  and  soil  is  moving  upward  at  ~40  cm/s  at  the 
surface . 

Figure  3.40  shows  the  air  overpressure  time  histories  at 
various  station  depths  for  Case  4.  Each  station  realizes  higher 
peak  air  overpressures  than  in  Case  3,  and  at  earlier  times.  The 
larger  particle  size  (and  thus  increased  permeability)  in  Case  4 
allows  the  air  to  permeate  the  soil  faster  and  to  greater  depths 
for  the  same  surface  pressure  loading.  This  allows  the  soil  to 
be  fluidized  at  earlier  times,  which  i3  seen  in  the  soil  velocity 
time  histories  in  Figure  3.41.  Equation  3.2  indicates  that  the 
earlier  time  leads  to  a  stronger  upward  pressure  force  at  that 
depth  for  the  fluidized  soil  in  Case  4.  The  resulting  upward 
soil  velocities  are  over  three  times  larger  than  calculated  in 
Case  3 . 


TIME  ISEC) 


Figure  3.41.  Soil  Velocity  Versus  Time  at  Various  Station  Depths 
for  Case  4  (D  =  0.142  cm). 


3.3  COMPARISON  OF  NUMERICAL  SIMULATIONS  WITH  EXPERIMENTS 


Although  the  DICE  simulations  were  not  performed  in  order  to 

3 

duplicate  the  laboratory  experiments  of  Zernow  et  al.  ,  qualita¬ 
tive  comparisons  can  be  made.  Both  are  one-dimensional  in  nature 
with  similar  pressure  loading  characteristics  at  the  soil/sand 
surface . 

Figure  3.2  shows  the  observed  surface  load  overpressure  time 
histories  used  for  the  two  types  of  laboratory  experiments  and 
the  overpressure  loading  function  used  in  Cases  1  to  4 .  While 
the  initial  peak  overpressure  is  6.89  bars  for  all  three,  the 
decay  in  time  varies.  The  exponential  decay  function  used  in  the 
numerical  simulation  has  a  decay  rate  which  puts  it  between  the 
1.36"  diameter  orifice  experiments  and  the  1.78"  experiments  up 
to  time  t  *  .600  s.  The  1.78"  experiments  have  the  fastest  decay 
rate  with  a  positive  phase  duration  of  -.450  s. 

In  the  laboratory  experiments,  the  bottom  of  the  sand  column 
(-6  ft)  was  restricted  from  moving  downward.  In  the  numerical 
simulation  a  continuum  boundary  condition  is  imposed  at  the  grid 
bottom  (-200  cm),  allowing  soil  and  air  to  flow  through  the 
bottom.  Case  1  showed  that  the  unfluidized  soil  essentially 
moves  downward  as  a  rigid  body  (see  Figure  3.17).  The  fluidized 
soil  can  then  be  thought  of  as  moving  upward  relative  to  the 
unfluidized  3oil. 

Figure  3.42a  shows  the  time  histories  of  the  positions  (or 
trace)  of  various  soil  marker  points  for  Case  1.  These  tracer 
points  move  with  the  local  soil  velocity  at  any  given  time.  Up 
to  t  *  .200  s  there  is  an  apparent  uniform  downward  translation 
of  the  tracer  points  throughout  the  column.  The  actual  soil 
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Figure  3.42*  Tracer  Point  Positions  Versus  Time  for  DICE  Case  1 

(D  =0.1  cm),  (a)  Unadjusted  and  (b)  Adjusted,  based 
on  -175  cm  Tracer  Position. 


surface  (represented  by  the  tracer  initially  at  0  cm  depth)  is 
depressed  7.5  cm  by  t  *=  .200  s.  The  depth  of  fluidization  is 
40  cm  at  this  time;  below  40  cm  all  tracer  points  have  moved 
8.4  cm  indicating  that  the  unfluidized  soil  behaves  uniformly  at 
the  modeled  depths. 

A  measure  of  the  fluidized  soil  motion  with  respect  to  the 
unfluidized  soil  can  be  obtained  by  measuring  the  tracer  point 
positions  against  the  tracer  point  which  started  at  175  cm  depth 
(and  remains  in  essentially  unfluidized  soil).  Figure  3.42b 
shows  the  tracer  positions  with  the  adjustment  made  to  all  tracer 
points.  The  depth  to  which  significant  motion  due  to  fluidiza¬ 
tion  occurs  is  now  seen  more  clearly.  The  surface  tracer  point 
moves  up  almost  25  cm  with  respect  to  the  unfluidized  soil;  the 
tracer  point  initially  at  100  cm  depth  moves  up  ~2  cm. 

Figure  3.43  shows  the  unadjusted  and  adjusted  surface  tracer 
point  positions  versus  time  for  Cases  1  to  4 .  The  same  adjust¬ 
ment  made  in  Case  1  is  also  made  in  Cases  2  and  3 .  For  Case  4 
there  was  no  adjustment  made  because  by  t  *  .400  s  the  entire 
modeled  soil  column  was  fluidized.  Figure  3.43  reveals  the 
essential  differences  between  the  DICE  cases.  Increased  soil 
permeability  (due  to  increased  characteristic  particle  size 
diameter)  results  in  larger  surface  displacements  due  to  the  pore 
air  phenomena.  In  the  absence  of  gravity  (Case  3)  there  is  also 
increased  surface  displacement  for  a  given  permeability. 

The  adjustment  of  the  tracer  positions  allows  comparison  of 
the  relative  motion  of  the  fluidized  soil  in  the  calculations  to 
the  motion  of  the  soil  observed  in  the  experiments.  Figure  3.44 
compares  the  surface  displacement  observed  in  the  experiments 
with  the  adjusted  surface  positions  calculated  in  Cases  1  and  2. 


and  (b)  Adjusted. 
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Figure  3.44.  Main  Surface  Displacement  Versus  Time  Observed  in 
Experiments  (Zernow  et  al.,  1973)  and  Calculated 
by  DICE  (Cases  1  and  2). 


Case  1  is  remarkably  similar  to  the  1.78"  test.  Both  surfaces 
start  to  move  upward  at  t  -.200  s  and  maintain  an  approximately 
constant  upward  velocity  until  t  -.400  s.  By  t  *  .450  s  the 
1.78"  tests  no  longer  have  any  surface  pressure  loading  (see 
Figure  3.2).  Since  there  is  no  longer  a  decrease  in  pressure  at 
the  surface,  the  upward  pressure  force  is  no  longer  sustained 
since  the  upward  moving  pore  air  serves  to  balance  out  the  pore 
air  pressure  differential.  Eventually  the  sand  motion  is 
dominated  by  gravity  as  the  air  can  no  longer  drag  the  fluidized 
sand  with  it.  By  t  -.550  s  the  sand  starts  to  fall  back.  On  the 
other  hand,  the  surface  pressure  loading  in  Case  1  is  still 
decaying  with  time  having  not  decayed  quite  as  rapidly  at  the 
1.78"  experiments  up  to  t  -.400  s.  Thus,  there  is  still  an 
upward  pressure  force  experienced  by  the  air  to  at  least 
t  -.600  s. 

The  surface  tracer  for  Case  2  also  started  to  rise  at 
t  -.200  s,  but  maintained  a  smaller  constant  upward  velocity  to 
t  - . 500  s.  However,  the  smaller  displacements  and  velocities 
compared  to  Case  1  are  indicative  of  a  smaller  sustained  upward 
pressure  force  resulting  from  the  permeability  dependence  of  the 
pore-air  phenomena. 

The  surface  displacements  for  the  1.36"  experiments  do  not 
start  to  rise  significantly  until  nearly  t  -.400  s  due  to  the 
slower  rate  of  decay  of  the  surface  pressure  loading  than  in  the 
1.78"  experiments.  However,  the  1.36"  experiments  also  exhibit  a 
constant  upward  velocity  for  a  substantial  length  of  time  due  to 
the  pore-air  phenomena. 

Figure  3.45  shows  the  position  of  the  near-surface  marked 
sand  layers  in  time  for  the  1.78"  experiments  compared  to  the 


Figure  3.45.  Comparison  of  Near-Surface  Marked  Sand  Layer  Positions  Versus  Time 

Observed  in  Experiments  (Zernow  et  al .,  1973)  and  Adjusted  Near-Surface 
Tracer  Point  Positions  for  DICE  Case  1  (D  =0.1  cm). 


equivalent  adjusted  soil  tracer  positions  for  Case  1.  Again, 
there  are  similarities  between  the  two  sets  of  curves.  There  is 
a  uniform  rise  of  the  near-surface  layers  beginning  after 
t  =  .200  s  and  continuing  beyond  t  =  .400  s.  In  the  1.78" 
experiments  the  sand  starts  to  move  upward  more  abruptly  and 
faster  than  in  Case  1,  probably  due  to  the  more  rapid  decay  rate 
of  the  surface  pressure  loading.  The  deeper  layers  for  both 
appear  to  start  moving  upward  at  slightly  later  times  than  those 
at  the  surface,  indicating  the  downward  progression  of  fluidiza¬ 
tion.  After  t  =  .500  s  the  soil  is  continuing  to  rise  for  Case 
1,  although  decelerating,  whereas  in  the  1.78"  tests  the  sand  is 
beginning  to  fall  downward  due  to  the  end  of  any  surface  pressure 
loading  at  t  -.450  s. 

Figure  3.46  compares  the  adjusted  near-surface  soil  tracer 
positions  of  Cases  1  and  2.  Evident  is  the  slower  rise  of  the 
soil  in  Case  2.  The  surface  tracers  begin  to  move  upward  just 
after  t  -.150  s  for  both  cases,  while  the  tracer  originally  near 
15  cm  depth  starts  moving  upward  near  t  ~.200  s  for  Case  1  and 
t  -.300  s  for  Case  2,  indicating  a  slower  downward  progression  of 
fluidization  in  Case  2. 

Figure  3.47  shows  the  time  of  fluidization  at  various  depths 
for  Cases  1  to  4 .  Fluidization  as  defined  here  occurs  when  the 
soil  particles  are  supported  entirely  by  air,  i.e.,  the  pressure 
associated  with  the  particles  supporting  each  other  is  zero.  For 
Case  1  the  "fluidization  front"  moves  down  through  the  soil  at 
~4  m/s.  The  fluidization  front  moves  more  slowly  at  m/s  for 
the  less  permeable  soil  modeled  in  Cases  2  and  3.  In  Case  3,  the 
soil  at  a  given  depth  is  fluidized  earlier  than  in  Case  2  because 
there  is  no  initial  lithostatic  pressure  distribution  in  the  soil 


Figure  3.46.  Adjusted  Near-Surface  Tracer  Point  Positions  Versus  Time  for  DICE 
Case  1  (D  =  0.1  cm)  and  Case  2  (D  =  0.05  cm). 
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to  be  overcome  by  the  air  before  fluidization  occurs.  The  highly 
permeable  Case  4  yields  a  fluidization  front  velocity  of  ~7  m/s. 

Some  measurements  of  overpressure  time  histories  were 
obtained  at  selected  depths  in  the  sand  column  for  the  laboratory 
experiments.  The  available  data  is  summarized  in  Table  3.1. 

These  can  be  compared  with  similar  data  from  Cases  1  to  4. 

Figure  3.48  shows  the  calculated  peak  air  overpressures 
versus  depth  for  Cases  1  to  4  compared  to  the  observed  experi¬ 
mental  data.  The  more  permeable  modeled  soils  (larger  particle 
size)  allow  more  air  to  penetrate  to  a  given  depth,  yielding 
higher  overpressures.  All  cases  display  an  exponential  decay 
with  depth  of  the  peak  overpressure  related  to  the  decay  of  the 
surface  loading.  The  less  permeable  modeled  soils  yield  more 
rapid  decay  rates  with  depth.  The  Case  2  and  Case  3  curves 
indicate  that  down  to  200  cm  depth  there  is  essentially  no 
gravity  dependence  on  the  peak  pore-air  overpressure  attained  at 
a  particular  depth. 

The  1.36"  standard  sand  tests  appear  to  exhibit  a  decay  in 
peak  overpressure  similar  to  that  of  Case  1;  the  observed  peak 
overpressures  are  slightly  lower  than  calculated  for  Case  1.  The 
observed  peak  overpressure  for  the  1.78"  tests  at  152  cm  depth  is 
comparatively  lower  than  for  the  1.36"  standard  sand  tests  at 
146  cm  depth.  With  the  more  rapid  decay  (and  thus  shorter 
duration)  of  the  surface  pressure  loading  for  the  1.78"  tests, 
less  air  is  forced  into  the  sand  column  than  for  the  longer 
duration  1.36"  tests  over  an  equivalent  length  of  time.  Less  air 
will  result  in  lower  peak  air  overpressures  at  a  given  depth  for 
sands  of  equivalent  permeability  and  porosity. 
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1  a  ted  (DICE  Cases  1  to  4)  and  Observed 
3)  Peak  Overpressure  at  Various  Depths 


For  the  same  surface  pressure  loading  the  peak  air  over¬ 
pressures  are  expected  to  be  larger  for  the  more  permeable  coarse 
sand  than  for  the  standard  sand.  At  the  146  cm  depth  this  is 
true.  The  83  cm  depth  data  point  for  coarse  sand  appears  to  be 

anomalous  in  this  regard,  and  is  left  unexplained  by 
3 

Zernow  ejfe.  al. 

Related  to  this  anomoly  could  be  the  drastically  reduced 
main  surface  displacement  for  the  coarse  sand  experiment  when 
compared  to  the  standard  sand  (see  Figure  3.4).  The  DICE 
calculations  assume  homogeneous  soil  mixtures  whose  perme¬ 
abilities  are  determined  by  the  soil  porosity  and  the  character¬ 
istic  particle  size  diameter.  The  calculations  showed  that 
increased  permeabilities  result  in  increased  surface  lofting  as 
well  as  higher  peak  pore-air  overpressures  at  a  given  depth.  In 
the  coarse  sand  experiments,  however,  a  breakdown  in  the  homo¬ 
geneous  3and  lattice  during  the  surface  pressure  loading  or 
during  subsequent  soil  expansion  could  cause  relatively  large 
diameter  air  channels  to  develop.  This  would  allow  the  pore-air 
to  vent  more  readily  and  directly  upward  through  the  channels. 
Thus,  there  would  be  a  decrease  in  the  sustained  upward  internal 
force  (drag)  available  to  loft  the  sand  upward. 

Figure  3.49  shows  the  calculated  peak  air  overpressure 
arrival  times  versus  depth  for  Cases  1  to  4  and  the  observed 
experimental  data  from  Table  3.1.  Figure  3.50  compares  the  times 
at  which  the  air  overpressure  first  achieves  10%  of  its  peak  at  a 
particular  depth.  The  more  permeable  modeled  soils  allow  air  to 
travel  faster  down  through  the  soil,  yielding  earlier  arrival 
times  for  both  the  peak  overpressure  and  10%  level.  The  slope  of 
the  peak  arrival  time  versus  depth  for  each  case  is  remarkably 
similar  to  those  in  Figure  3.47  which  shows  the  soil  fluidization 


123 


4— 

O 

-X 

CD 

03 

E 

CD 

•r- 

CL  — * 

i— 

o 

o 

r-.  lo  o  i 

4-  CD 

LO  to  OO  1 

i — 

o  to 

f-H 

^■<*1  •  tmmMmt: 

03 

_ _ <• 

• 

> 

•r— 

o 

S- 

r-H 

03 

03  o  < — 
or  ro 

oo 

a)  cn 

CD 

E  4-  ' 

E 

•r-  CD 

•f— 

h-  4->  * 

h- _ 

CD  . 

</)  E 

03 

11 

o 

CO  OsJ  CNJ  ( 

3  IO  tj 

CD 

CO 

C\J  CM  CO  < 

</)  *r- 

a. 

• 

•  •  • 

^.0  4^ 

•  P— 

CD  0) 

S- 

>  = 

S- 

00  3 

0)  o 

s-  s-  c; 

c 

c/1  4-  CD 
t/1  M 

CD  (/) - 

s-  x: 

CD 

II 

CL4->  (/> 

i- 

S-  Q.-t-> 

3 

<u  <u  C 

o> 

>Q  CD 

CO  to 

ID  CO  04  i 

O  E 

03 

CD  s- 

CO 

co  rs  co  < 

TD  -r- 

CD 

5-  « 

• 

«  •  • 

4-  CD  $- 

a. 

Q-JD 

*“H  »-H 

O  P  d) 

t-  — 

U  Q- 

CD 

>*>  CD  X 

> 

S-  i—  UJ 
ro  CD 

O 

E  00  CD 

p  .o 

§  +->  Z3 

00  03  J— 

x: _ 

a3  co  ro 

• 

c= 

• 

•  •  • 

r—l 

a-  g 

C\J 

co  co  to  i 

• 

cu  J-L 

LO 

CO  CO  1 

co 

Q 

1—1 

CD 

O 

03 

h— 

CD 

^ v 

a. 

-O 

“D 

>> 

CD  C 

CD  C 

O  03 

O  03 

•r-  00 

•r-  OO 

-*-> 

4- 

4-  1 

c 

•r-  “O 

•r-  *U 

CD 

s-  s- 

S-  i- 

E 

O  03 

O  03  1 

•r“ 

“O 

T3 

=  c 

=  C 

CD 

00  03 

CO  03 

Q. 

*-> 

CO  4-> 

X 

•  OO 

•  OO 

UJ 

r-H 

1.36"  Orifice  83.8  1.03  .20  .056 

(Coarse  Sand)  146.3  .81  .275  .100 


CALlFflRNlA  RESEARCH  AND  7ECHNSL8GT.  INC. 
a  ICE  3280 

TIME  OF  ARRIVAL  OF  UK  PEAK  OVERPRESSURE  VS.  DEPTH 


Calculated  (DICE  Cases  1  to  4)  and  Observed 
,  1973)  Time  of  Arrival  of  10%  Level  of  Peak 


time  versus  depth.  For  Case  1  the  peak  air  overpressure  appears 
to  move  down  through  the  soil  at  ~4  m/s.  The  less  permeable 
soils  of  Cases  2  and  3  yield  ~2  m/s,  while  the  highly  permeable 
Case  4  yields  -6  to  7  m/s.  For  each  case  there  is  also  a  delay 
of  -.080  to  .120  s  after  the  arrival  of  the  peak  overpressure  at 
a  given  depth  until  soil  fluidization  begins  at  that  depth.  The 
1.36"  experiments  also  yield  earlier  arrival  times  for  the  more 
permeable  coarse  sand  as  compared  to  the  standard  sand.  The  data 
fall  around  the  curves  for  Case  1  but  the  slopes  of  lines  through 
the  data  would  indicate  faster  velocities  for  the  peak  over¬ 
pressure  and  10%  level. 


SECTION  4 


MISERS  BLUFF  11-2 

4.1  EXPERIMENTAL  RESULTS 

The  MISERS  BLUFF  1 1-2  experiment  involved  the  simultaneous 
detonation  of  six  high-explosive  charges  (120  ton)  forming  a 
hexagon  100  m  on  a  side.  Figure  4.1  shows  the  experimental 
configuration  and  the  measurement  gauge  locations  for  data  used 
in  this  report.  Details  of  the  experiment  can  be  found  in 
Reference  1. 

Figure  4.2  shows  the  overpressure  time  history  (after 
detonation)  measured  at  various  azimuthal  angles,  all  25  m  from 
the  center  of  the  hexagon.  Three  separate  peaks  are  evident  at 
each  angle,  although  not  as  distinct  at  0°.  The  time  histories 
for  the  two  gauges  at  30°  and  210°  are  very  similar.  By 
t  -  .150  s  the  positive  phase  has  ended  for  each,  and  a  small 
negative  phase  is  sustained  until  t  ~.400  s. 

Figure  4.3  shows  various  vertical  velocity  time  histories  of 
ground  motion  gauges  originally  at  0.5  m  depth  and  25  m  radius. 
The  azimuthal  angles  are  different  from  those  for  the  overpres¬ 
sure  time  histories.  Figure  4.4  shows  the  corresponding  vertical 
displacement  of  each  gauge  as  a  function  of  time.  The  velocities 
and  displacements  are  very  similar  for  each  gauge.  Up  to 
t  -.125  8  the  ground  motion  is  downward,  corresponding  to  the 
airblast  pressure  load  pushing  down  at  the  surface.  Peak 
downward  velocities  are  5  to  10  m/s.  From  t  -.125  s  to  t  -.450  s 
the  ground  motion  is  upward,  with  peak  velocities  of  ~5  m/s  at 
t  -.300  s.  Peak  upward  displacements  of  ~1  m  are  achieved.  This 
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0° 


270°. 


Figure  4.1.  Location  of  Charges  and  Various  Measurement 
Gauges  in  MISERS  BLUFF  1 1-2. 
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Figure  4.3.  Vertical  Ve  .ocities  Measured  in  MISERS  BLUFF  1 1 -2 

from  Ground  Motion  Gauges  Originally  at  0.5  m  Depth. 
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Figure  4.4.  Vertical  Displacements  Measured  in  MISERS  BLUFF  1 1 -2 
from  Ground  Motion  Gauges  Originally  at  0.5  m  Depth. 
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sustained  upward  motion  is  similar  to  the  motion  seen  in  the 
laboratory  experiments  described  in  Section  3;  however,  now  the 
upward  velocities  are  associated  with  negative  pressures.  From 
t  -.450  s  to  t  -.850  s  the  ground  motion  is  again  downward,  after 
which  time  the  velocities  oscillate  about  zero. 

4.2  NUMERICAL  SIMULATION 


A  numerical  calculation  (Case  5)  was  performed  with  the 
multiphase  DICE  code  in  order  to  examine  the  apparent  pore-air 
phenomenology  seen  in  the  ground  motion  data  of  MISERS  BLUFF 
1 1 -2.  The  calculation  was  designed  using  one-dimensional  (1-D) 
planar  geometry  in  the  vertical  direction,  with  a  surface 
overpressure  loading  function  defined  by  an  approximation  to  the 
airblast  measured  at  25  m  radius  from  the  center  of  the  array. 

Figure  4.5  shows  the  overpressure  loading  function  used  in 
Case  5,  along  with  the  airblast  measured  at  30°  and  210°  for 
comparison.  There  are  three  distinct  peaks  which  occur  in  the 
loading  function  before  t  -.110  s  after  detonation;  the  magni¬ 
tudes  are  6.5  bars  at  t  *  .045  s,  12.8  bars  at  t  -  .063  s,  and 
7.7  bars  at  t  -  .102  s.  The  measurements  showed  larger  peaks 
occurring  a  couple  of  milliseconds  earlier  than  the  assumed 
loading  function.  A  negative  phase  occurs  after  t  -.150  s  with 
overpressures  reaching  -0.8  bars.  By  t  -.400  s  the  overpressure 
loading  function  becomes  zero. 

Figure  4.6  shows  the  recommended  calculat ional  soil  profiles 
from  Reference  11  developed  to  characterize  the  subsurface 
layering  beneath  the  120°  and  300°  radials  of  MISERS  BLUFF  1 1-2. 
Figure  4.7  shows  the  same  for  the  single  charge  event  MISERS 
BLUFF  II-l,  held  ~400  m  away  from  event  2.  It  is  this  profile 
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Figure  4.5.  Overpressure  Loading  Function  Versus  Time  Used  in 
DICE  Case  5  and  Airblast  Time  Histories  Observed 
in  MISERS  BLUFF  11-2. 
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which  was  used  in  determining  the  soil  layer  thicknesses  used  in 
Case  5. 

Figure  4.8  shows  the  1-D  grid  used  in  the  calculation. 

Three  soil  layers  are  used  to  model  the  site.  The  top  soil  layer 
consists  of  five  0.12  m  cells.  Each  subsequent  cell  in  Layer  2 
is  10%  larger  than  the  cell  above  it,  down  to  7.35  m  depth. 
Thereafter  the  cells  are  uniform  (0.73  m)  down  through  Layer  3. 

Table  4.1  defines  the  MISERS  BLUFF  1 1-2  material  properties 
used  to  model  each  of  the  three  soil  layers.  Each  soil  layer  is 
modeled  using  only  one  characteristic  particle  size,  which  is 
estimated  from  data  found  in  Reference  12.  The  permeability 
model  discussed  in  Section  2.3  yields  permeability  as  a  function 
of  porosity  and  particle  size.  The  remaining  material  properties 
are  from  Reference  11.  The  initial  densities  increase  with 
depth,  with  the  corresponding  initial  porosities  decreasing  with 
depth  from  44%  to  1%.  In  addition,  hysteresis  effects  are 
included  and  the  modeling  is  described  on  Table  4.1.  The  loading 
and  unloading  constrained  moduli  for  each  layer  were  determined 
from  the  stress-strain  curves  in  Reference  11. 

The  initial  modeled  permeabilities  for  Case  5  are  much 
smaller  than  in  Cases  1  to  4.  In  Cases  1  to  4  the  initial 
permeabilities  were  in  the  lO2  -  103  Darcy  range.  In  Case  5,  the 
initial  permeabilities  are  different  in  each  modeled  soil  layer  - 
k—.04  Darcy  for  the  top  layer,  k  ~70  Darcy  for  Layer  2  and 
k~.003  Darcy  for  Layer  3.  One  of  the  results  of  a  decrease  in 
permeability  was  seen  in  the  Case  3  and  Case  4  comparison,  which 
demonstrated  that  the  air  would  permeate  the  soil  slower  and  to 
shallower  depths  for  the  same  surface  pressure  loading  and 
equivalent  porosities.  In  effect,  the  internal  interaction 
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Soil  Layer  1 


(Dp  =  .0005  cm) 


;7D 


Soil  Layer  2 
(Dp  =  .04  cm) 


Soil  Layer  3 
(Dp  =  .08  cm) 


Table  4.1.  MISERS  BLUFF  II-2  Material  Properties. 


Vertical  Strain 


forces  between  the  air  and  soil  particles  (i.e.,  drag)  are 
increased  with  a  decrease  in  permeability,  so  that  the  relative 
velocities  will  be  smaller. 

Case  5  includes  the  effects  of  a  uniform  gravitational  field 
(g  =  980  cm/s2).  This  requires  an  initial  lithostatic  soil 
pressure  at  depth  to  support  the  weight  of  the  total  amount  of 
soil  above  that  depth.  The  initial  pore-air  pressure  was  taken 
to  be  1  atm,  uniform  in  each  layer. 

Figure  4.9  to  4.13  shows  the  pressure  and  velocity  distri¬ 
butions  for  soil  and  air  from  t  =  .047  s  after  detonation  to 
t  =  .600  s.  The  first  airblast  peak  (6.5  bars  overpressure)  in 
the  loading  function  occurs  at  t  =  .043  s.  Figure  4.9  shows  that 
the  relative  velocity  between  the  soil  and  air  is  essentially 
zero,  due  to  the  low  permeability  of  the  top  soil  layer.  This 
contrasts  with  the  large  relative  velocities  seen  in  Case  1  (see 
Figure  3.9)  shortly  after  the  initial  application  of  the  surface 
pressure  load.  For  Case  5,  the  soil  is  supporting  virtually  all 
of  the  overpressure  load,  except  for  a  small  amount  supported  by 
the  air  which  has  barely  begun  to  permeate  the  soil  near  the 
surface . 

Figure  4.10  shows  that  by  t  *  .101  s  all  modeled  soil  layers 
have  been  affected  by  the  airblast,  which  is  at  its  third  and 
final  peak  of  the  positive  phase.  The  soil  and  air  are  moving 
downward  together  -  almost  1  m/s  at  the  surface  and  ~0.4  m/s  at 
13  m  depth.  At  the  interfaces  between  each  soil  layer  (~0.6  m 
and  ~4.8  m  depths)  discontinuities  in  the  material  properties 
cause  propagation  of  reflected/transmitted  pressure  and  velocity 
waves . 
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DEPTH  IN)  46.7b  MSEC 

(a)  Soil,  Air  and  Total  Stress  Versus  Depth,  and 

(b)  Soil  and  Air  Velocity  Versus  Depth  at 
t  =  .047  Seconds  for  Case  5  ( MB  1 1 - 2 ) . 


Figure  4.9. 
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DEPTH  tM)  0.203  SEC 

Figure  4.11.  (a)  Soil,  Air  and  Total  Stress  Versus  Depth,  and 

(b)  Soil  and  Air  Velocity  Versus  Depth  at 
t  =  .203  Seconds  for  Case  5  (MBII-2). 


PRFSSURF  (BARS) 
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Figure  4.12.  (a)  Soil,  Air  and  Total  Stress  Versus  Depth,  and 

(b)  Soil  and  Air  Velocity  Versus  Depth  at 
t  =  .400  Seconds  for  Case  5  ( MB  1 1 -  2 ) . 


By  t  =  .203  s  the  airblast  is  in  its  negative  phase  (-0.8 
bar  overpressure).  Figure  4.11  shows  that  all  of  the  soil  in 
Layer  1  and  part  of  Layer  2  (to  ~2.0  m  depth)  is  f luid ized .  The 
soil  and  air  are  moving  upward  above  ~1  m  depth  in  the  soil  with 
a  peak  velocity  of  over  12  m/s  at  the  surface;  the  surface  has 
risen  ~U .  5  m  by  this  time. 

For  DICE  Cases  1  to  4,  fluidization  at  depth  occurred  as  a 
result  of  enough  air  permeating  the  soil  for  the  pore-air  to 
support  the  positive  (decaying)  overpressure  load  applied  at  the 
surface.  The  same  effect  is  achieved  when  a  negative  overpres¬ 
sure  load  is  applied,  even  in  the  absence  of  a  positive  phase. 

The  negative  overpressure  permits  the  pore-air  to  support  the 
reduced  pressure  load.  And,  the  pore-air  expands  to  reach 
pressure  equilibrium  with  the  applied  surface  "underpressure". 

By  t  =  .400  s  the  surface  overpressure  loading  has  returned 
to  zero.  Figure  4.12  shows  an  air  pressure  minimum  at  ~1  m  above 
the  original  ground  surface.  The  resulting  pressure  gradient 
causes  downward  acceleration  for  the  near-surface  material. 

Peak  downward  velocities  of  almost  5  m/s  exist  at  the  surface, 
which  has  been  raised  almost  2  m  above  the  original  surface. 

Peak  upward  velocities  of  ~7  m/s  are  present  at  -1.0  m  above  the 
original  ground  surface.  The  difference  in  velocities,  although 
still  small,  between  the  soil  and  air  in  the  0  to  2  m  depth 
region  occurs  in  soil  Layer  2,  where  the  permeability  is  10 3 
times  larger  than  in  Layer  1. 

The  surface  overpressure  load  remains  at  zero  after 
t  »  .400  s.  Figure  4.13  shows  that  by  t  =  .600  s  the  pressure 
minimum  has  moved  down  to  ~1.5  m  depth,  where  the  peak  upward 
velocities  are  ~2  m/s  for  the  soil  and  air.  A  peak  downward 


velocity  of  over  4  m/s  is  seen  at  the  surface.  All  of  soil 
Layer  2  is  moving  upward  at  this  time;  below  the  interface 
between  Layer  2  and  3  the  velocities  are  relativey  small. 

The  use  of  an  approximate  airblast  function  based  on  data 
measured  in  MISERS  BLUFF  1 1-2  permits  direct  comparison  with 
measured  around  motion  data.  In  Figures  4.14  and  4.15  the 
velocity  and  displacement  time  histories  of  the  tracer  point 
originally  at  0.5  m  depth  are  compared  with  the  ground  motion 
data  measured  with  the  four  gauges  originally  at  0.5  m  depth. 
Phenomenologically  the  agreement  appears  good.  Up  to  t  '-.150  s 
the  calculation  and  data  show  downward  motion,  though  the 
calculated  motion  is  not  as  large  as  measured.  This  could  he  due 
to  the  smaller  peak  pressures  assumed  in  the  loading  function 
(see  Figures  4.2  and  4.5).  Both  curves  show  a  sustained  upward 
velocity  (in  the  negative  overpressure  phase)  from  t  ~ . 150  s  to 
t  ~.400  s,  with  the  calculated  velocities  and  displacements  being 
somewhat  larger  than  measured.  After  this  time  the  motion  is 
again  downward  for  both.  Results  from  the  Case  5  calculation 
indicate  the  primary  role  of  pore-air  phenomena  in  causing  the 
sustained  upward  ground  motion  data  measured  in  MISERS  BLUFF 
1 1  -2 , 
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01CE  3280  N1SERS  BLUFF  11-2  R-2S  M 

HEASURED  ANO  CRLCULA7E0  TRACER  POINT  VELOCITY  VS.  TIKE  ZO-O.S  H 


TIME  ISEC) 

Vertical  Velocity  Time  History  Comparisons,  Case  5  Versus 
MISERS  BLUFF  1 1 -2  Experimental  Data  at  25  m  Radius  and  0.5 
Initial  Depth. 
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ATTN:  FCTT ,  J.  Lee 
ATTN:  FCTT 
ATTN:  FCTXE 
ATTN:  FCTT,  W.  Summa 


Joint  Chiefs  of  Staff 

ATTN:  GD10,  0-5  Nuc 


&  Chem  Div 


DEPARTMENT  OF  THE  ARMY  (Continued! 

US  Army  Engineer  Ctr  l  Ft  Bel  voir 
ATTN:  Tech  Library 

US  Armv  Engrg  Div,  Huntsville 
ATTN:  HNDEN-FO 

US  Army  Engrg  Waterways  Exper  Station 
ATTN:  WESSC,  J.  Jackson 
ATTN:  J.  Zelasko 
ATTN:  Library 
ATTN:  WESSA,  W.  Flathau 
ATTN:  WESSS,  J.  Ballard 
ATTN:  WESSE,  D.  Day 

US  Army  Material  &  Mechanics  Rsch  Ctr 
ATTN:  Tech  Library 

US  Army  Materiel  Dev  &  Readiness  Cmd 
ATTN :  DRXAM-TL 

US  Army  Nuc  i  Chem  Agency 
ATTN:  Library 
ATTN:  MONA-WE 
ATTN:  MONA-WE,  J.  Uecke 


DEPARTMENT  OF  THE  NAVY 

Naval  Civil  Engrg  Lab 

attu.  i  AT.  J.  Forrest 


Naval  Postgraduate  School 

ATTN:  Code  1424  Library 
ATTN:  G.  Lindsay 


Joint  Strat  Tgt  Planning  Staff 
ATTN:  JPPFM 
ATTN:  JPSS 
ATTN:  JPTM 
ATTN :  JLKC 
ATTN:  JLKS 
ATTN:  JLK,  DNA  Rep 

Under  Secy  of  Def  for  Rsch  &  Engrg 

ATTN:  Strat  *  Thtr  Nuc  Force,  B.  Stephan 

2  cy  ATTN:  Strat  &  Space  Sys (OS) 

DEPARTMENT  OF  THE  ARMY 

BMD  Advanced  Technology  Ctr 
ATTN:  ATC-T 

BMD  Systems  Command 

ATTN:  BMDSC-HLE,  R.  Webb 


Chief  of  Engineers 

ATTN:  OAEN-MPE-T 

Harry  Diamond  Labs 

ATTN:  DELHD-NW-P 
ATTN:  DELHD-TA-L 

US  Army  Ballistic  Rsch  Labs 

ATTN:  DRDAR-BLT ,  0.  Keefer 
ATTN :  DROAR-BLA-S 

US  Army  Cold  Region  Res  Engrg  Lab 
ATTN:  Library 


Naval  Rsch  Lab  ,  ,  .  ,  . 

ATTN:  Code  2627,  Tech  Lib 
ATTN:  Code  4040,  J.  Boris 
ATTN:  Code  4040,  D.  Book 

Naval  Surface  Weapons  Ctr 
ATTN:  Code  F31 
ATTN:  Code  X211,  Tech  Lib 
ATTN:  Code  R44,  H.  Glaz 


Naval  Surface  Weapons  Center 

ATTN:  Tech  Library  &  Info  Svcs 


Br 


DEPARTMENT  OF  THE  AIR  FORCE 

Air  Force  Institute  of  Technology 
ATTN:  Library 
ATTN :  ENA 


Air  Force  Systems  Command 
ATTN :  DLWM 


Air  Force  Weapons  Lab 

ATTN:  NTE,  M.  Plamondon 

ATTN:  NTES-A 

ATTN :  SUL 

ATTN:  NTED 

ATTN:  NTES-C 

ATTN:  NTED-S 

ATTN:  NTES-R 

ATTN:  NTEO 

ATTN:  DEY 


DEPARTMENT  OF  THE  AIR  FORCE  (Continued) 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 


Air  University  Library 
ATTN:  AUL-LSE 

Assistant  Chief  of  Staff,  Intel  1 
ATTN :  IN 

Assistant  Chief  of  Staff 
Studies  &  Analysis 
2  cy  ATTN :  AF/SAMI ,  Tech  Info  Div 

Asst  Secy  of  the  AF 
Rsch,  Dev  A  Logistics 

ATTN:  SAFALR/Oep  for  Strat  &  Space  Sys 

Ballistic  Missile  Office/DAA 


ATTN: 

EN 

ATTN: 

ENSN 

ATTN: 

EN,  C 

Case 

ATTN: 

ENSN, 

A.  Schenker 

ATTN: 

ENSN, 

E.  Furbee 

Deputy  Chief  of  Staff 
Rsch,  Dev  A  Acq 

ATTN:  AFRO-M,  Spec  Asst  for  MX 
ATTN :  AFROQI 

Strategic  Air  Command 


Applied  Rsch  Associates,  Inc 
ATTN:  R.  Frank 

Applied  Rsch  Associates,  Inc 
ATTN:  0.  Piepenburg 

Applied  Rsch  Associates,  Inc 
ATTN:  J.  Shinn 

Applied  Theory,  Inc 
2  cy  ATTN:  J.  Trulio 

ARTEC  Associates,  Inc 
ATTN:  S.  Gill 

Astron  Rsch  A  Engrg 

ATTN:  J.  Huntington 

Boeing  Co 

ATTN:  Aerospace  Library 
ATTN:  S.  Strack 

California  Rsch  A  Technology,  Inc 


ATTN 
ATTN 
2  cy  ATTN 
2  cy  ATTN 


Library 
K.  Kreyenhagen 
M.  Rosenblatt 


ATTN : 

DEPM 

2  cy  ATTN 

D. 

Orphal 

ATTN: 

NRI/STINFO  Library 

ATTN: 

XPFS 

California 

Rsch 

t  Tech,  Inc 

ATTN: 

XPQ 

ATTN 

F. 

Sauer 

OTHER  GOVERNMENT  AGENCY 

Central  Intelligence  Agency 
ATTN:  OSUR/NED 

DEPARTMENT  OF  ENERGY  CONTRACTORS 

University  of  California 

Lawrence  Livermore  National  Lab 
ATTN:  D.  Glenn 

Los  Alamos  National  Laboratory 
ATTN:  M.  Sandfor 
ATTN:  C.  Keller 

Sandia  National  Laboratories 

ATTN:  Org  7112,  A.  Chabai 
ATTN:  Div  1111,  J.  Reed 
ATTN:  0.  Bannister 

DEPARTMENT  OF  DEFENSE  CONTRACTORS 

Acurex  Corp 

ATTN:  C.  Wolf 

Aerospace  Corp 

ATTN:  H.  Mirels 

ATTN:  Library  Acq  Ml/ 199 

Agbabian  Associates 

ATTN:  M.  Agbabian 

Applied  Rsch  Associates,  Inc 
ATTN:  N.  Higgins 
ATTN:  J.  Bratton 


Carpenter  Rsch  Corp 

ATTN:  H.  Carpenter 

University  of  Denver 

ATTN:  Sec  Officer  for  J.  Wisotski 

Dynamics  Technology,  Inc 
ATTN:  D.  Hove 

Electro-Mech  Systems,  Inc 
ATTN:  R.  Shunk 

H&H  Consultants,  Inc 
ATTN:  W.  Hall 
ATTN:  J.  Haltiwanger 

H-Tech  Labs,  Inc 

ATTN:  B.  Hartenbaum 

I IT  Research  Institute 

ATTN:  Document  Library 

Karnan  AviDyne 

ATTN:  R.  Ruetenik 

Karnan  Sciences  Corp 

ATTN:  E.  Conrad 

Karnan  Tempo 

ATTN:  DAS I AC 

Karnan  Tempo 

ATTN:  DAS  I AC 

Martin  Marietta  Denver  Aerospace 
ATTN:  D-6074,  G.  Freyer 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 


DEPARTMENT  OF  DEFENSE  CONTRACTORS  (Continued) 


Kendan  Associates 

ATTN:  R.  Dannenberg 
ATTN :  J .  Owren 

Maxwell  Laboratories,  Inc 
ATTN :  J .  Murphy 

McDonnell  Douglas  Corp 
ATTN:  H.  Herdman 
ATTN:  D.  Dean 
ATTN:  R.  Halprin 


Rand  Corp 

ATTN:  B.  Bennett 

Rand  Corp 

ATTN:  P.  Davis 

S-CUBED 

ATTN:  C.  Needham 


Science  &  Engrg  Associates,  Inc 
ATTN:  H.  Linnerud 


Merritt  CASES,  Inc 

ATTN:  Library 


Science  &  Engrg  Associates,  Inc 
ATTN:  B.  Chambers  III 


Mission  Rsch  Corp 

ATTN :  C .  Longmi re 

University  of  New  Mexico 
ATTN:  G.  Leigh 
ATTN:  J.  Kovama 
ATTN:  P.  Lodde 
2  cy  ATTN:  D.  Calhoun 

Pacific-Sierra  Rsch  Corp 
ATTN:  D.  Gormley 

Pacific-Sierra  Rsch  Corp 

ATTN:  H.  Brode,  Chairman  SAGE 


Science  Applications,  Inc 
ATTN:  W.  Layson 
ATTN:  J.  Cockayne 
ATTN:  G.  Binninger 

Science  Applications,  Inc 
ATTN:  H.  Wilson 
ATTN:  Tech  Library 
ATTN:  R.  Schlaug 

SRI  International 

ATTN :  0 .  Col  ton 
ATTN:  G.  Abrahamson 
ATTN:  Library 


Pacifica  Technology 

ATTN:  Tech  Library 
ATTN:  R.  Allen 

Patel  Enterprises,  Inc 
ATTN:  M.  Patel 

Physical  Rsch,  Inc 

ATTN:  R.  Deliberis 


Physics  International  Co 
ATTN:  Tech  Library 
ATTN:  R.  Collins 

R&D  Associates 

ATTN:  C.  Knowles 
ATTN :  J .  Lewi s 
ATTN:  Tech  Info  Ctr 
ATTN:  A.  Kuhl 
ATTN:  D.  Simons 
ATTN :  P .  Haas 

S-CUBED 

ATTN:  C.  Dismukes 
ATTN:  K.  Pyatt 
ATTN:  Library 
ATTN:  J.  Barthel 


Teledyne  Brown  Engineering 
ATTN:  F.  Leopard 
ATTN:  B.  Hartway 
ATTN:  D.  Ormond 

Terra  Tek,  Inc 

ATTN:  Library 
ATTN :  J .  Schatz 

TRW  Electronics  A  Defense  Sector 
ATTN:  N.  Lipner 
ATTN:  Tech  Info  Ctr 

TRW  Electronics  &  Defense  Sector 
ATTN:  G.  Hulcher 
ATTN:  E.  Wong 
ATTN :  P .  Dai 

Weidlinger  Associates,  Consulting  Engrg 
ATTN:  T.  Deevy 

Weidlinger  Associates,  Consulting  Engrg 
ATTN:  I.  Sandler 

Weidlinger  Associates,  Consulting  Engrg 
ATTN:  J.  Isenberg 

R&D  Associates 

ATTN:  G.  Ganong 
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DTIC 


